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ABSTRACT 


The Li K and Mg L,, spectra have been observed from a range of Li-Mg 
alloys, evidence for electron transfer from Mg being obtained. More distinct 
changes have been observed for the same spectra when Li was evaporated 
onto solid Mg and Mg onto solid Li. Changes in the Li K and Al L,, were 
observed when Li was evaporated onto Al and in the Al L,, when Al was 
evaporated onto Mg. No changes greater than 5 x 10-2ev were observed in 
edge wavelengths or edge widths although marked depopulation of higher 
energy states occurred in Mg and Al. It is inferred that changes in the screening 
of K and L levels:compensate almost exactly for the effects of electron 
transfer. The results of Catterall and Trotter (1959) were confirmed in that 
the Mg band contracted and the Li band extended on the low-energy side. 
Li deposited on Al shows a K spectrum without any metal edge, suggesting 
that the first zone is filled and that there is no overlap with the second zone. 


§ 1. INTRODUCTION 


THE effect of alloying on the valence band energy levels of metals of low 
atomic number has been the subject of a number of investigations. Earlier 
work seeking possible changes in the soft x-ray spectra was hampered 
by the difficulties of photographic photometry in this region which makes 
accurate measurements of intensities almost impossible to obtain. The 
attendant need to dissipate large amounts of power in the alloy target 
and the inability of the photographic method to detect the short-term 
effects of target contamination are added handicaps. As a consequence 
few definite and confirmed effects resulting from alloying have been 
recorded. 

Photoelectric detection based on photon-counting enables the repro- 
duction and measurement of intensities to 1 or 2% in this region. The 
power dissipated in the target is reduced to about one-hundredth of that 
needed for photography, and changes of intensity with time resulting 
from target contamination can be observed and therefore often avoided. 

The first changes caused by alloying in the Li-Mg system observed by 
photoelectric means were reported by Catterall and Trotter (1959). 
Evidence was obtained pointing to electron transfer from the constituent 
of higher to the constituent of lower valency while the spectrum of each 
constituent retained a recognizable though somewhat modified identity. 
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The observations reported below, though based initially on the same 
alloy system, recorded the most marked changes in metallic combinations 
not normally regarded as involving alloying since they consisted of 
evaporated layers of one constituent on a solid, cooled base of the other. 
Since the changes can be most effectively interpreted in terms of electron 
transfer, alloying, even though it be very restricted in the dimension of 
depth, is assumed to occur. 


§ 2, EXPERIMENTAL 


The emission spectra were obtained with the photon-counting spectro- 
meter previously described (Fisher et al. 1958). In the initial experiments 
on Li-Mg alloys, target currents up to 11ma and voltages of 4kv were 
used. The pressure in the specimen chamber was about 10-°mm and 
carbonaceous contamination of the target surface by decomposed residual 
vapours was serious. With this type of contamination intensities decrease 
linearly with time (Crisp 1958a). Subsequently the vacuum in the speci- 
men chamber was improved to 1-2 x 10-*mm and the target was almost 
completely surrounded by a copper baffle cooled to —115°c by external 
conduction (Pollard 1953) to a rod immersed in liquid air. The cooling 
could be maintained indefinitely and the target was kept under vacuum 
for one or two days before making observations. So long as the baffle 
was cooled below —100°c carbonaceous contamination was reduced to a 
negligible amount. A glass window allowed observation of the interior 
of the specimen chamber. 

The target current was stabilized to 0-5% (Crisp 1958b), enabling 
relative intensities to be reproduced to within 1 or 2%. The variable 
speed drive was replaced by a gearbox with fixed speeds so that spectra 
could be accurately reproduced to a chosen wavelength scale. The use of 
more sensitive photomultipliers} allowed the reduction of target currents 
toabout l!ma. The aluminium grating used in previous work was replaced 
by a glass grating of similar constants. The resolution, calculated accord- 
ing to the criterion of Fisher (1954), was maintained at 0-9 A. 

A four-sided target with faces of Cu, Al, Mg and Be was used so that 
spectra for the relevant elements either from the solid or from evaporated 
layers on Cu could be recorded without breaking the vacuum. Solid Li 
specimens were accommodated in a different target. The tungsten carbide 
scraper was used to clean target surfaces when required. The target 
area covered by the electron beam is about 10x 5 mm. 

All spectra reproduced have been traversed from the high-energy side 
so that the recorder trace (i.e. the time axis) runs from right to left. The 
chart speed was always one inch per minute. If contamination is occurring 
the low-energy side will be reduced in intensity and the maximum may 
be shifted slightly towards higher energies. All wavelengths in a given 
spectrum appear to show similar rates of contamination with the exception 
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of the Li ‘line’ at 2534 which grows approximately as [1 —exp(—a2t)]. 
The appearance of this ‘line’ and the closely associated fluroescence of the 
Li target will be discussed elsewhere, but it should be recorded that after 
it was discovered that a massive evaporation of Li has the effect of suppress- 
ing the appearance of the ‘line’ this precaution was always taken before 
recording a Li spectrum. 

Our LiK spectrum has been published elsewhere (Crisp and Williams 
1960). The MgL,, and AlL,, spectra seen by our instrument are shown in 
figs. 2 and 5} respectively. Although the resolution is insufficient to record 
the full height of the Al ‘spike’, the small kink at about 1814 is always 
seen and constitutes a real feature of the spectrum not previously recorded. 
Its position corresponds to the minimum in Matyas’ theoretical N(E) 
curve (‘Tomboulian 1957) which marks the lowest energy for states in the 
second Brillouin zone. 

Edge wavelengths, taken at half maximum intensity of the edge were 
checked by reference to the fiducial marks made on the recorder trace at 
every revolution of the lead screw. Variations less than 0-2 4 (equivalent 
to about 0-05ev) are considered to be within the limits of experimental 
error. Edge widths, taken between the points of 5°4 and 95% maximum 
intensity can be measured to the same accuracy. 


§ 3. Liratum—MAGNeEsium ALLOYS 


Solid specimens of alloys of Mg—Li containing 15, 40 and 70at.% Li 
were examined. The alloys were received as metal bars contained in a 
clear, but unidentified packing oil. Specimens were milled to shape using 
the packing oil as cutting fluid. They were pressed into a frame which was 
soldered to the target holder with Woods’ metal. Heating of the specimens 
was in this way kept to a minimum. Any oxidation of the target surface 
was subsequently removed in vacuum by the scraper. Thermal contact 
of the specimens with the water-cooled target holder was poor. 

At target currents up to 11 ma and a pressure 8 x 10-mm a rapid time 
variation of intensity was observed for both Mg and Li. The form of the 
variation is shown in fig. 1. As the Li intensity increased a strong blue 
fluorescence appeared over the bombarded area. ‘The time variation 
suggested that under the action of the beam Li atoms were migrating to 
the surface. Although spectra can be recorded at a time when the variation 
of intensity is slow, it is doubtful whether the surface layer examined has 
the same composition as the original alloy. 

Subsequently, in the improved vacuum and with improved sensitivity, 
the intensity variation was examined for target currents between 0-7 
and 4ma. The intensity was then found to decay exponentially for both 
constituents with a half-intensity period about one minute for the 30/70 
and 60/40 alloys and somewhat larger for the 85/15 alloy. This is approxi- 
mately the same as the half-period for growth of the Li ‘line’. 


SE EE 2 re See eee 
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Figure 2 shows recorded spectrat for the three alloys together with the 
spectrum of pure Mg. Contamination in these cases is slower but in ine 
30/70 alloy the Mg has disappeared in less than 9 min and has been replace 

i ‘line’. 
ereenee the changes in the Mg spectra agree with the results of Be 
and Trotter (1959). There is a distinct diminution of the Mg spike ; 
indicating a depopulation of the higher electron states. This diminution 
becomes more marked as the fraction of Li increases, though its exact 
extent is masked by the overlap of the Li ‘line’. The edge wavelength 
does not change appreciably and there is no convincing evidence that the 


Fig. 1 
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Time variation of intensity (counts per second) for Li K and Mg L,, spectra in 
various Li-Mg alloys. Target current 11 ma voltage 4kv. Pressure 
10-> mm. 


edge is broadened. The low-energy maximum appears to shift from 2614 
in pure Mg to about 257A in the alloys, but this may again be caused by 
the Li ‘line’. Because of low intensities and the effects of contamination 
it was not possible to obtain information about the shape of the Li K 
spectrum. 


§ 4. EvaporatEep LAYERS 


Since it appeared, from the intensity changes obtained with high target 
powers, that Li atoms were migrating toward the surface and possibly 
producing a lithium-rich layer on the surface, it was decided to look at the 
Mg spectrum when a lithium-rich layer was deliberately produced by 
evaporating Li onto a solid Mg target. The result was surprising. The 
Mg spectrum was modified to a much greater extent than in the alloy 


SUE EERE ee tr re ters f 

{ Date, target, full-scale count rate, percentage error at half-full scale (Elmore 
and Sands 1949) target current and voltage, optical resolution and pressure 
are recorded on each spectrum reproduced. 
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spectra. Subsequently, Mg was evaporated onto Li, producing further 
big changes in the Mg spectrum. The effect of evaporating Li, Na and K 
onto Be, Mg and Al was examined. Combinations for which alloying 
does not normally occur showed no spectrum changes. The only appreci- 
able changes involved combinations of Li, Mg and Al which are discussed 
below. 

Examples of recorded spectra illustrating the effect of evaporating Li 
onto Mg and Al, Mg onto Li and Al onto Mg are shown in figs. 3, 4, and 5. 
Contamination was appreciable only in the case of Mg on Li where it has 
the effect of reducing intensities on the low energy side of the spectrum. 

The changes in the spectra are so marked that it is assumed they are due 
to alloying, either instantaneously by diffusion following release of latent 
heat by the condensing element, or very rapidly under the action of the 
bombarding electrons, even though the underlying constituent is held 
at the temperature of the cooling water. In the interval of less than a 
minute between ceasing evaporation and recording the maximum of the 
Mg peak the changes have been completed. Alloying in thin films has 
been discussed by Weaver and Hill (1959) who record that Michel (1956) 
observed immediate alloying following successive evaporation of Ag—Zn 
and Ag-Sn. Michel’s suggestion that alloy formation is observable for 
metal pairs whose diffusion constants at room temperature are larger 
than 10-*° cm?/sec is not consistent with the observation of changes in the 
L,3 spectrum of Al deposited on Mg, the only pair for which diffusion 
constants are available. Fujiki (1959) has observed instantaneous 
alloying following successive depositions of metal vapours on a Formvar 
film, but in this case the latent heat released is sufficient to raise the tempera- 
ture of the film. However, Fujiki showed that alloy formation was 
promoted by small grain size and large area of the substrate surface. Target 
surfaces subjected to repeated scraping would meet these requirements. 
The Li-Mg, Li—Al and Al-Mg systems have not been studied by electron 
diffraction. 

Film thicknesses could not be quantitatively controlled, although a 
degree of control was obtained by running the furnace at minimum power 
for evaporation and switching off at different intervals after a visible film 
first appeared. In this way ‘thin’ and ‘thick’ films could be deposited. 
As the emission spectrum is produced in a surface layer of the order of 
104A thick (Skinner 1940), it is assumed that for appreciable emission 
to be observed the film thickness could be at least of the order of 10° A. 

The spectra of each element are discussed separately below to facilitate 
comparison of the modifications with the spectrum of the pure metal. 


§ 5. MAGNESIUM 


The L,3 spectrum of pure Mg seen by our instrument is shown in fig. 2. 
In fig. 3 are shown the spectra of Mg when overlain by Li and when 
evaporated onto solid Li. In fig. 6 various Mg spectra have been super- 
imposed with the intensity at the low-energy maximum adjusted to a 
constant value, the remainder of the spectrum being normalized accordingly. 
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Three spectra of pure Mg were averaged ; three of Mg under Li, four of Mg 
over Li and four of Mg in Mg,Sn. The wavelength of the high-energy edge 
has in each case been calculated by comparison with the edge of pure Mg 
obtained with the grating in the same adjustment. The averaged curves 
are taken from spectra in which evaporation has been continued until no 
further changes occur in the shape of the Mg spectrum. A light evaporation 
reduces the Mg spike slightly making the spectrum resemble that from an 
Mg-Li alloy. 


Fig. 6 
> 
Magnesium L543 Spectra 
sas eee praeen~, | f 
le ta aca 
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Smoothed averages of three spectra of pure Mg, three of Mg under Li, four of Mg 
deposited on Li and four of Mg in powdered Mg,Sn. Intensities nor- 
malized to make the low-energy maxima equal. 


From fig. 6 it is concluded that the edge wavelength is constant within 
the experimental uncertainty (5x 10-ev), and that the width of the 
edge between the 5% and 95% maximum intensity points is similarly 
constant. Although the change in position of the low energy maximum 
from 260A for pure Mg to about 2624 for Mg under Li is not established 
with certainty, the change to a flatter maximum with a mean at 2674 
for Mg on Li is real. This change represents a decrease of more than one 
electron volt in the energy of the states making contact with the boundary 
of the first Brillouin zone. Together with this shift in the maximum, the 
ratio of maximum intensity at the edge to maximum intensity in the band 
is 68% (probable error 1-5% from 17 spectra) for Lion Mg and 5 1% (probable 
error 3% from four spectra) for Mg on Li. The two spectra are quite 
distinct. On the low-energy side the Mg bandwidth is decreased in Mg 
under Li relative to pure Mg by about 0-5ev. A similar decrease is found 
by comparing averaged spectra taken with the aluminium grating. This 
is in agreement with the results of Catterall and Trotter for Mg-Li alloys. 
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Although the Mg on Li shows a more marked decrease, this is largely 
if not wholly due to contamination. Assuming the contamination produced 
an exponential decrease in intensity (since it was associated with the 
appearance of fluorescence and the Li ‘line’) a rough correction of the 
intensities in the tail shows no definite decrease in the band width. 

Li can dissolve in Mg to about 18 at. % and since the Mg remains in its 
h.c.p. state the Mg spectrum should not be much changed in shape. This 
is found to be so for Lion Mg. Mg can dissolve in Li to about 70 at. % Mg 
but the structure is b.c.c. and the Mg spectrum should be modified 
appropriately. The distinct difference between the Mg spectra noted 
above suggests that the Mg deposited on Li dissolves in the Li. 

Between the maximum at 2624 and the kink at 2554 the intensities 
for pure Mg are higher than for Mg under Li. This difference, together 
with the possible shift of the maximum may result from a dilution of the 
spectrum of Mg containing dissolved Li by a small contribution from the 
spectrum of Mg dissolved in Li. In the case of Mg on Li the situation is 
reversed and the spectrum of Mg on Li may be diluted though the distinctive 
shape of the maximum is not lost. 

Figure 6 also includes an averaged Mg spectrum from the intermetallic 
compound Mg,Sn in powder form. Counting rates were low and the 
background was high in the four spectra averaged and only the general 
trend of the curve is significant, the position of the low energy maximum 
and the breadth of the band near the tail being to some extent dependent 
on the choice of background intensity. Since Mg,Sn has a fluorite structure 
it was to be expected that the Mg L,, spectrum from this compound would 
be modified. However, it closely resembles the spectrum of pure Mg, the 
edge wavelength and edge width being unchanged. There is some indication 
of an increased band-width on the low-energy side. The result does not 
support Farineau’s (1937) observation of similar K spectra for the con- 
stituents of the Al-Mg intermetallics, Al,Mg, and Al,;Mg,. 

When Al is evaporated onto Mg no changes are observed in the Mg 
spectrum although as will be seen from fig. 5 the AlL,, is changed. Since 
a normal Mg L,, spectrum will be produced from the solid Mg target, this 
will be observed as well as that from the Mg in contact with the Al. 


§ 6. ALUMINIUM 


Figure 4 shows three recordings of the AlL,, spectrum for increasing 
thicknesses of Li deposited on solid Al. Figure 5 shows the normal Al 
spectrum seen by our instrument in which the resolution is insufficient 
to show the full extent of the high energy spike. There is also shown in 
fig. 5 the spectrum of Al evaporated onto Mg. Figure 7 shows the result 
of averaging three spectra of Al, four spectra of Al with a superimposed 
‘thick’ layer of Li and two spectra of Al deposited on Mg. 

The point chosen on the Al spectrum at which to normalize the intensities 
has been taken at about 1894. This is well below the maximum intensity 
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of states in the first zone and was chosen on the assumption that the 
population of states in the first zone of Al would not be affected by alloying 
because the second zone is half-filled. This assumption is not borne out 
when the results are analysed. 

The trend of absorption by Li is also indicated in fig. 7 following Skinner 
and Johnson (1937), with an arbitrary zero corresponding to the absorption 
at the low-energy side of the K edge. Li absorption should be roughly 
constant to the high-energy side of the normalizing wavelength. On the 
low-energy side there are two increases to a maximum at 1964 which 


Fig. 7 


Aluminium Lo3 Spectra 
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Smoothed averages of three L,, spectra of pure Al, four spectra of Al with a 
superimposed ‘ thick ’ layer of Li and two spectra of Al deposited on Mg. 
Intensities have been normalized to the same value at 189 A. Kinks in 
the pure Al spectrum at 173 A and 181 Aare reproducible. Li absorption 
in this region is also shown with an arbitrary zero and scale. 


qualitatively correspond to the two regions in which the intensity of the 
modified spectrum of Al falls below that of pure Al. However, no quan- 
titative match has been achieved and it appears that the N(E) curve of Al 
is modified in this region. The spectrum of Al on Mg shows distinct 
divergencies from that of pure Al in the same region and since these latter 
cannot be attributed to absorption they constitute further evidence for 
changes in the N(E) curve for the first zone. 

The prominent maximum at 183A may have been enhanced by the 
choice of normalizing wavelength and by Li absorption. Its position 
corresponds to the maximum density of states in the first Brillouin zone 
of pure Al. To the high-energy side the shape of the spectrum resembles 
Matyas’ (Tomboulian 1957) theoretical N(E) curve with excision of 
states of higher energies. If, for the ‘thick’ layer of Li, it is assumed that 
no further change in the Al spectrum will result from increasing the 
thickness of Li, then by analogy with the case of Li on Mg, no dilution by 
the spectrum of normal Al from the underlying solid target should be 
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present. Therefore the modified Al spectrum is the result of a combination 
of (i) structure changes which produce a different N(E) distribution, (ij) 
depopulation of higher energy states of the first and second zones resulting 
from electron transfer, and (iii) Li absorption whose undetermined 
extent makes estimates of the decrease in band width at the low energy 
side very hazardous. 

For Al on Mg a normal Al spectrum produced in the unalloyed outer 
region of the Al film will dilute that from modified Al. The Al spike has 
not disappeared as with the thick Li layer and there is no absorption. 
However, structure changes in the alloyed Al could produce changes in 
the N(E) curve at lower energies. The choice of 189 A as the normalizing 
wavelength then loses its significance. 

The edge wavelength and edge width are apparently not changed and 
no edge shift to 171-8A reported by Gale and Trotter (1956) for Mg—Al 
alloys is present. Nothing can be said about the band width. There is 
no resemblance in the Al on Mg spectrum to the alloy spectra of Gale 
and Trotter. 

§ 7. Liraium 


The K spectrum of Li shows a distinct high energy edge extending to 
about 75% of the intensity at the maximum (Crisp and Williams 1960). 
On the low-energy side the shape and width of the band can be seriously 
affected by development of the Li ‘line’ whose maximum comes at 253A 
(Fisher e¢ al. 1958). After the discovery that the development of the 
line and the associated fluorescence could be stopped by massive evapora- 
tion of Li as a ‘getter’ this was always done before recording Li spectra. 
The Li spectra of figs. 3 and 4 and the smoothed curves of fig. 8 are therefore 
not affected by the Li line. 

Figure 8 shows the result of averaging four spectra of pure Li, four 
spectra of Li on Mg and four spectra of Li on Al. These have been 
normalized to the same intensity at the maximum and for the first two 
the edge is taken to extend to 75% of the maximum. Only two spectra of 
Li under Mg were recorded, neither suitable for analysis, and all that can 
be said is that for Li under Mg no differences from normal Li have been 
established. 

Figure 3 shows that when Li is deposited on Mg there is a relative 
increase in intensities on the low-energy side of the spectrum. This 
accords with the observations of Catterall and Trotter for Li-Mg alloys. 
There is no observable change in edge wavelength, edge width or energy of 
the maximum. The latter energy is determined by the contact of the 
anisotropic Fermi surface with the first zone boundary and that this 
should be unchanged suggests that emission from the overlying b.c.c. Li 
is observed as well as that from the Li presumably dissolved in the h.c.p. 
Mg. Catterall and Trotter obtained a similar result for an alloy containing 
20-9% atomic Li which presumably is mainly a solution of Li in Mg. 

Figures 4 and 8 show that when Li is deposited on Al the Li spectrum 
is changed very radically. The high-energy edge has disappeared, the 
maximum has shifted and the band width has increased. The ‘edge 
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width’ to 75% maximum intensity is now about 1-2ev, comparable with 
the ‘edge width’ shown by non-conductors (Skinner 1940). ‘The maximum 
has shifted about 0-5ev to lower energies. The band width, including 
the tail, is about 8ev compared with 4-2ev for pure Li, for which the 
extrapolated band width is 3-0 + 0-1ev (Crisp and Williams 1960). If this 
spectrum comes from Li atoms at the interface combined with Al as 
LiAl, the first Brillouin zone of Li could be practically full. A spectrum 
such as that seen would be expected if the zone was full and an energy 
gap existed between first and second zones. Any zone overlap would 
prevent the disappearance of the edge. 


Fig. 8 


Lithium K Spectra 
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Smoothed averages of four K spectra of pure Li, four spectra of Li on Mg and 
four spectra of Li on Al. Intensities normalized to make the maxima 
equal. Dotted curves show a hypothetical separation of the Li on Al 
into two components, assuming the low-energy bulge is symmetrical 
about the maximum at 247 A. 


It should be emphasized that the low-energy bulge has nothing to do 
with the Li ‘line’. Its maximum is at 2474 compared with 2534 for 
the ‘line’. No fluorescence, which always accompanies the line, was 
visible on the target when these spectra were obtained. Figure 8 shows a 
tentative separation of the Li spectrum into two components assuming 
that the low-energy bulge is symmetrical about its maximum. The low 
energy component might result from the appearance of excitation states 


(Parratt 1959) of Liin LiAl. The width of the high energy component is 
about 5ev. 


§ 8. Discusston 


The spectrum changes noted in Li-Mg alloys are seen to a much greater 
extent with evaporated layers. This is taken as evidence for instantaneous 
alloying at the interface. Observations made by Catterall and Trotter 
(1959) on Li-Mg alloys are supported by our results for both alloys and 
evaporated films. There is no change greater than 5x 10-2ev in edge 
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wavelengths, or edge widths. The band width for the constituent of 
higher valency is lessened, the change showing up on the low-energy side. 
The band width of the constituent of lower valency is increased, the 
increase showing on the low-energy side. The exception is the case of 
Li deposited on Al in which the band shape is radically changed. 

Since the Mg L,, under Li shows a complete loss of the spike, which is 
the result of second zone overlap, it is assumed that the depopulation of 
energy states due to electron transfer has progressed at least as far as the 
energy at the base of the spike, emptying the second zone entirely and the 
first zone also to some extent. However, the absence of any measurable 
change in energy of the edge shows that this emptying of the upper levels 
has been offset very effectively by changes in the screening of the L levels. 
The screening compensation is so much more exact than for the case of 
copper—zine alloys discussed by Varley (1956), that a close connection 
between the extent of electron transfer, leading to changes in the popula- 
tion of energy states and the changes in screening of the inner levels, is 
indicated. 

Since the edge has, in effect, moved up, the tail follows, and the energy 
range over which depopulation has occurred should be estimated from 
the narrowing of the band at the tail. It will be seen from fig. 6 that this 
energy range is approximately equal to that from the edge to the foot of 
the spike. If the structure remains the same the low-energy maximum 
should shift to higher energies by the same amount. This is not observed. 
As pointed out above, it is possible that the spectrum of h.c.p. Mg contain- 
ing dissolved Li is diluted by the spectrum of Mg dissolved in b.c.c. Li 
for which the maximum is shifted to lower energies by the structure 
change. The required shift of low-energy maximum is seen in fig. 3 of 
Catterall and Trotter (1959) but the maximum for their pure Mg occurs at 
264A compared with 260 +14 in our fig. 6. 

The above considerations cannot apply to the intermetallic Mg,Sn 
because it has a fluorite structure and the N(E) curve should be different. 
However, the increased band width, if real, is seen on the low-energy side 
as in Li, the edge energy being unchanged. 

As regards Lion Mg, there is agreement with the observations of Catterall 
and Trotter that the intensities increase in the tail while the maximum 
does not changeenergy. Ifshielding is reponsible for there being no change 
in edge energy, increased population of high-energy states should move 
the maximum to lower energies and decrease the height of the edge. This 
has not been observed, either by Catterall and Trotter, for an alloy contain- 
ing 20-9% atomic Li or by ourselves for Lion Mg. Whilst in our case the 
Li spectrum from the unalloyed part of the Li film would dilute that from 
the modified Li, the lack of any change at all in the maximum makes this 
an inadequate explanation. 

The case of Li deposited on Al involves the most drastic change recorded 
in a soft x-ray spectrum. The energy at half-intensity of the modified 
Li spectrum is very nearly equal to the edge energy in pure Li, and the 
shift in maximum is approximately the same as the low-energy increase 
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in band width. This accords with the arguments stated above. However, 
the shape of the maximum has changed considerably although the structure 
change from Li to LiAl is slight. If the tails and the maxima of the two 
spectra are brought into coincidence the pure Li edge approximately 
bisects the modified Liband. This supports the inference that the complete 
loss of high energy edge indicates a completely filled first zone with an 
energy gap above it. 

The case of Al is more difficult to interpret. If, as assumed above the 
evaporation of Li was carried so far that no dilution by the normal AlL,3 
spectrum from the underlying solid is present, the observed Al spectrum 
results from the three effects previously listed. A structure change from 
f.c.c. Al to b.c.c. LiAl would produce a modification of the N(E) curve for 
Al, but its true extent is disguised by Li absorption. 

Until the spectrum from bulk LiAl can be observed it is futile to analyse 
the modified Al spectrum. Dilution from pure Al may be reponsible for 
the greater prominence of the high energy spike in the Al overlying Mg. 
However, if this is so the modified Al spectrum must possess a more 
prominent maximum at 1824 and the modification of the low-energy 
region would be more severe than appears in fig. 7. 
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ABSTRACT 


Thin foils of 99-999% aluminium prepared by electropolishing have 
been found to develop surface pits along the grain boundaries during ageing 
at room temperature. It is suggested that the pits are formed by a vacancy 
mechanism, and that they are associated with the initiation of etch attack 
along the grain boundaries of aged pure aluminium. 


§ 1. INTRODUCTION 


TuE localization of etching attack of pure aluminium has been widely 
associated with the emergence of dislocations (see for example Amelinckx 
1953, Lacombe 1947, Lacombe and Beaugard 1948) or with the segregation 
of impurities to the dislocation lines (Wyon and Crussard 1951, Forty 
and Frank 1955). In the case of the etching of slip lines on strained 
aluminium, Bassett and Edeleanu (1960) have advanced an alternative 
suggestion, that it is the rupture of the natural oxide film and its sub- 
sequent failure to repair perfectly, which plays an important part in 
the initiation of etch attack. 

Braun, ef al. (1958) have found that periods of room temperature 
ageing are necessary after the preparation of the specimen surface, 
before grain boundary etching is observed. They have suggested this 
behaviour is due to diffusion of impurities along the grain boundaries 
and into the oxide film. 

A consideration of the likely residual impurities in high purity 
aluminium, and the thermodynamics of the replacement of aluminium 
in the oxide film with these impurities, makes this view difficult to 
accept. The work described now suggests a possible alternative 
explanation of grain boundary attack after ageing treatments. 


§ 2. EXPERIMENTAL 
Specimens of polycrystalline aluminium of 99:999% purity were 
prepared as foils thin enough for transmission electron microscopy, 
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by electropolishing in a bath of one part perchloric acid to five parts 
ethyl alcohol. The initial foil was 0-001 in. thick and had a highly 
polished surface finish. Prior to electropolishing foils were heated to 
600°c followed by air cooling. Electron microscope examination was 
carried out using a Siemens Elmiskop I operated at 80 or 100 kv, 
specimens being mounted in the stereo holder so that dislocations in 
the foils could be brought into satisfactory contrast. 

After thinning, foils were aged, mounted on grids ready for examination 
in the microscope, in a vacuum desiccator for periods up to four months. 
Some foils were given an accelerated ageing treatment by maintaining 
at a temperature of 100°c in an air oven for 24 hours. For comparison 
with the thin film observations anodic film replicas of the foil surfaces 
have also been made. 


§ 3. RESULTS AND INTERPRETATION 


Some typical grain boundaries in a foil aged at room temperature 
for four months after thinning are shown in fig. 1 (a) and (b)7, and at a 
higher magnification in fig. 2. These micrographs reveal the formation 
of small surface pits along the grain boundaries. This pit formation 
is quite general, all grain boundaries appearing pitted. The higher 
magnification micrograph shows one example in which the pits are 
bounded by well-formed crystallographic faces. For comparison speci- 
mens before ageing have been examined and found to have grain 
boundaries which are completely free of pits of this type. 

Accelerating the ageing process by holding a thin foil at 100°c for 
12—24 hours it has been possible to photograph in the electron microscope 
the same grain boundary both before and after ageing. Figures 3 and 4 
show such a pair of micrographs in which the appearance of pits at the 
grain boundary at A and B is clearly seen. A piece of dirt on the left- 
hand side of the picture forms a convenient point of reference. 

The size of pits which have been observed lies between 1004 and 
10004; some tendency for the pits to be larger but fewer in specimens 
aged at room temperature has been noticed compared with those in 
which the ageing was accelerated. Figure 5 gives some information 
of the depth of the pits where an array can be seen at both foil surfaces 
on an inclined grain boundary. The depth of these pits is about one-fifth 
of the foil thickness, judged from the interception of the extinction 
contours in the grain boundary. This figure makes quite evident that 
the pits are confined to the foil surfaces. 

The relative shallowness of the pits has made their detection on oxide 
replicas uncertain. Although some oxide replicas have shown evidence 
of the pitting along the boundaries this has not always been seen. 
Carbon replicas have been tried but in general have failed to detect the 
pits. An explanation could be that the pits have formed under the 
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natural oxide film and a carbon replica would then not reveal the pits 
unless the film was completely collapsed. 

The electron micrograph, fig. 6, is typical of an aged thin foil which 
has been very lightly etched in HCl/HNO,/HF. The disposition of 
the holes made in the foil by this treatment strongly suggests that etching 
started at the grain boundary pits produced on ageing. This type of 
etching behaviour was confined to aged specimens and was not observed 
on freshly prepared ones. 


§ 4. Discussion 


Thin film electron microscope observations have shown that after 
ageing surface pits form along the grain boundaries in aluminium. Such 
behaviour can best be interpreted as a migration and aggregation of 
vacancies at the surface of metal and beneath the natural oxide film. 
The formation of surface pits on aluminium has previously been 
reported by Doherty and Davies (1959). In their experiments the pits 
visible under an optical microscope formed after a small, but critical, 
drop of temperature, after holding the aluminium specimen at a high 
temperature. Doherty and Davies interpret their results in terms of 
the aggregation of the excess vacancies produced by the temperature 
drop. It is interesting to note that they observe that pits form only in 
the grains and that near grain boundaries no pits are found; the bound- 
aries acting as very effective sinks for vacancies. 

In the present experiments any possible concentration of thermal 
vacancies quenched in from the heat treatment is too small to account 
for the pits. It is probable in any case that many quenched-in vacancies 
are annealed out during electropolishing when a thin foil may reach 
quite high temperatures. The present pits presumably must also result 
from segregation of vacancies but it would seem more reasonable that 
in this case the vacancies migrate across the specimen surface or grain 
boundary and that metal atoms move away from the region of the pit. 
The driving force is presumably the large surface energy associated 
with the T-like junction made by the grain boundary and outer surface 
and is the same as in thermal etching. The final stable condition in 
that case would be a trench rather than pits as now found but parametric 
or crystallographic factors may favour the formation of pits. It is not 
possible from the present work to say whether the effect is associated 
with a grain boundary diffusion of vacancies and metal or whether the 
outer surface is involved, but it seems likely in either case that the film 
above the voids will be disturbed and it is suggested that this will then 


favour local etching. 
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ABSTRACT 


Etching of ice crystals grown from both the vapour and the melt has 
revealed etch pits in concentrations of up to 108/em? which apparently mark 
the sites of emergent dislocations. Silvered Formvar replicas show the etch 
pits on the basal surface to be hexagonal pyramids on which spiral or closed 
terraces of steps are often visible. These features are interpreted in terms of 
the etching of helical dislocations. The etch pits tend to group themselves 
along low-angle grain boundaries and slip lines. When the crystal is 
mechanically stressed, etch channels appear on the basal plane along the 
< 1120) directions which suggests that these are the preferred glide directions. 

Double rows of small etch pits radiating from impurity specks and local 
centres of high stress probably mark the sites of dislocation loops inter- 
secting the surface. Pyramidal hillocks appearing in concentrations of 
10°/em? on the basal faces of crystals grown from the vapour may be growth 
features and not the result of etching. 


§ 1. LyrRopUCcTION 


THERMAL or chemical etching acts by the preferential evaporation or 
dissolution of the surface at imperfections where molecules are surrounded 
by fewer nearest neighbours than usual and so are less strongly bound. 
The shape of an etch pit will depend upon the geometrical arrangement 
of the molecules in the particular crystal face and upon the relative rates 
of dissolution in directions parallel and perpendicular to the surface. 
The geometry of the etch pits may therefore allow the crystal face and the 
directions of crystallographic axes to be identified. 

But much more important, it has recently been demonstrated, by several 
workers, that etch pits mark the locations of emergent dislocations. This 
has been most convincingly demonstrated by Jones and Mitchell (1957) 
who were able to decorate and make visible dislocations in silver-halide 
crystals by the preferential deposition there of metallic silver, and show 
that the terminations of dislocations in the crystal surface coincided with 
the positions of etch pits. Further evidence of this close association was 
provided by Johnston ef al. (1957) who were able to match the etch pits 
appearing on opposite sides of an etched thin slice of a silicon crystal, the 
inference being that the members of a matched pair were joined by a 
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dislocation line passing through the crystal normal to the surfaces. Further- 
more, Johnston and Gilman (1959), during successive stressing and etching 
of lithium fluoride crystals, observed that successive rows of etch pits 
formed in the manner expected if they originated on dislocations which 
moved progressively under the applied stresses. 

The rather extensive studies of etched crystals, mainly germanium and 
silicon, have revealed that : 


(i) Although it is by no means certain that all emergent dislocations 
are revealed by the appearance of etch pits (dislocations of large Burgers 
vector will obviously be favoured), observations of the pits provide useful 
information on the density and distribution of dislocations. 


(ii) Although screw and edge dislocations produce etch pits of similar 
shape, those associated with screw dislocations are generally smaller, more 
shallow, and less symmetrical than those associated with dislocations of 
strong edge character. 


(iii) Having revealed the existence of helical dislocations by the decora- 
tion of fluorite crystals, Amelinckx, Bontinck and Dekeyser (1957) have 
suggested that etching of such dislocations may explain the occasional 
appearance of spiral pits on germanium and silicon. 


(iv) Etch pits tend to group themsleves along slip lines and low-angle 
grain boundaries. The presence of etch pits not only reveals these features 
but a detailed study of the distribution, density, asymmetry and depth 
of the pits enables one to deduce the angle between grain boundaries and 
the directions of slip for both edge and screw dislocations. 


Because a study of etch-pit formation appeared to offer a simple and 
powerful method of studying the dislocation structure of crystals, we 
decided to exploit it in our studies of ice crystals where dislocations have 
not been revealed by the appearance of spiral terraces of growth steps nor 
by attempts to decorate them. 

The appearance of etch pits on the surfaces of ice crystals was first reported 
by Truby (1955). Hexagonal pyramidal features of base diameter }-20 pu 
and depth $-3 were observed and sometimes exhibited a stepped-layer 
structure concentric with the c axis. Higuchi (1958) and Muguruma and 
Higuchi (1959) discovered that etch pits could be produced on the surfaces 
of ice and snow crystals by etching with a Formvar solution. Hexagonal 
pits appeared on the basal surface of a snow crystal in concentrations of 
up to 10°/em®. After prolonged etching the bottom surfaces of the pits 
revealed the layer-type structure noted by Truby, the height of individual 
steps varying between 3 and 16. These Japanese authors have used the 
etch pits to determine the size distribution of, and orientation of the crystal 
axes in, individual grains of polycrystalline ice. 

The object of the present investigation was to make a more detailed 


study of etch pits on ice crystals and to relate them to the dislocation 
structure. 
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§2. Ercuine TECHNIQUE 


Surfaces of ice crystals, grown from both the vapour and the melt, 
were etched with a 1% solution of polyvinyl formal (Formvar) in ethylene 
dichloride. A small quantity of the solution was poured onto the ice 
surface and the solvent allowed to evaporate and leave a thin, transparent, 
plastic film on the surface which could be peeled off and examined as a 
replica. Etching occurred because traces of water in the ethylene dich- 
loride froze to produce filamentary ice crystals which subsequently 
evaporated to leave small holes in the Formvar film. Etching, by evapora- 
tion, then proceeded beneath these holes to produce wide, shallow pits 
on the ice surface. Smaller, pyramidal etch pits subsequently appeared 
inside these large pits as shown in fig. lj. Subsequent applications of 
the Formvar solution gave rise, in similar manner, to further large pits 
which merged with earlier ones so that, eventually, the whole surface 
became etched. The etch pits could be studied at any stage by peeling 
off the plastic film with a piece of transparent adhesive tape, fixing it 
face upwards on a microscope slide and viewing it under the microscope 
in transmitted light.. To obtain better contrast and resolution a thin 
layer of silver was evaporated onto the replica which, when viewed in 
reflected light with magnifications of x 400, revealed the detailed structure 
of the etch pits as shown in figs. 4 and 6. 


§3. Ercu Pits on Icz CrystaLs GROWN FROM THE MELT 


A rectangular metal channel was filled with distilled water and cooled 
slowly below 0°c in a cold room to produce a block of ice whose flat surfaces 
contained ice crystals of various orientations which could be identified 
by the shapes of the large etch pits. The ice block was placed in a 
refrigerator held at —20°c, its surface polished with fine emery paper 
followed by a silk cloth and plastic replicas of its surface taken at intervals 
of 24 hours. 

The first replica showed only the fine scratches produced during polishing 
but the second replica, taken 24 hours later, showed a flat surface covered 
by wide, shallow etch pits about 100 across with a surface density of 
about 103/em?. Etch pits appeared on all the exposed crystal faces, 
their shapes revealing the crystalline symmetry and the orientation of the 
face. Figure 2 shows shallow hexagonal pits on the basal surface, each 
about 100, in diameter, viewed in transmitted light. 

The replica in fig. 1 shows many small pyramidal etch pits (seen inverted), 
5 to 10 in diameter, formed inside a large, shallow pit on a basal surface. 
The smaller pits, appearing in concentrations of up to 10°/cm?, were 
hexagonal in cross section with their sides usually parallel, but occasionally 
perpendicular, to the sides of the large pit. The heights of the pyramids 
were roughly equal to the diameters of their bases. Under vertical illumina- 
tion the sides of the smaller pits appeared smooth but, when illuminated 
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at a slight angle to the vertical, they appeared to be built up of several 
concentric steps each a few tenths of a micron in height. On many of the 
larger pyramids (d > 10), these were accompanied by one, and sometimes 
two, large concentric steps of height up to ly. A few per cent of these pits 
exhibited hexagonal spirals of similar step height—see fig. 3. 

On some of the later replicas from completely etched surfaces, areas of 
basal surfaces were crossed by lines of hexagonal pyramidal pits forming 
etched channels. <A silvered replica showing an (inverted) line of pyramidal 
pits is shown in fig. 4. These were formed at a grain boundary, the average- 
spacing of the pits divided into the intermolecular distance giving 30” for 
the angle of misorientation. 

Ice crystals were also grown from the melt by slowly cooling a wide 
shallow dish of water in such a manner that freezing started at the surface 
to produce large, flat, single crystals. These were removed with tweezers, 
frozen to a microscope slide, and examined by transmitted light under a 
microscope in a cold room held at —8°c. Etching of the surface with 
successive applications of Formvar solution was carried out as before 
but now evaporation proceeded much more rapidly and the etch pits 
on the basal surface, both large and small, tended to be circular rather 
than hexagonal. However, concentric rings about ly in height and 
occasional spirals with one or two turns were still observed although smaller 
steps were now not visible. 

In an attempt to study the effect of the inclusion of small, insoluble, 
foreign particles in the ice, the shallow dish of water was allowed to collect 
dust for a few days. Ice crystals were then formed on the surface as 
before and silvered replicas of their basal surfaces, taken at 30 min intervals, 
now revealed many patterns similar to those shown in fig. 5. Radiating 
from a small pyramidal etch pit are six double etch channels which, having 
obtained a length of about 10, broaden to about double their original 
width and finally terminate in either another pyramidal etch pit or in 
the side of the large, shallow pit which contains them all. Broadening of 
the channel is sometimes accompanied by a change in direction of about 
30°. 


§4. SURFACE STRUCTURE OF CRYSTALS GROWN FROM THE 
VAPOUR 


Single ice crystals were grown from the vapour, under carefully controlled 
conditions of temperature and supersaturation, in the diffusion cloud 
chamber described by Hallett and Mason (1958). The crystals were 
grown on a fine nylon fibre hanging vertically down the centre of the 
chamber; the crystal habit varied along the length of the fibre from 
hexagonal plates>needles—prismatic columns->plates+>dendrites—>plates 
prisms as the temperature varied from 0°c to —50°c. Hexagonal plates 
of up to 5mm diameter, needles up to 2 cm long x 2mm wide, and dendrites 
3cm x lem were grown in this way and transferred to the cold room for 
etching and surface examination by the techniques described above. 
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Replicas of previously unetched basal faces of plate-like crystals 
revealed, instead of etch pits, raised hexagonal pyramids, usually less 
than 10,,in diameter. Under high magnification these exhibited concentric 
rings and spirals similar to those observed in etch pits. These hillocks, 
shown in fig. 6, were eroded by subsequent etching. They were accompanied 
by numerous smaller, circular imperfections too shallow to be definitely 
identified as either hills or depressions. Their surface density is about 
10°/cm?. 

Replicas of previously unetched prism faces showed a high density, 
about 10°/em?, of features which were rectangular, aligned parallel to 
the c axis, and had four sloping walls—see fig. 7. Again they were so 
shallow that it was difficult to determine whether they were hills or 
depressions. 

Etching with Formvar solution produced hexagonal etch pits on the 
basal surfaces of plate-like ice crystals in concentrations of order 104/cm?, 
the edges of the pits being parallel to the crystallographic (1120) direc- 
tions—see fig. 8. Pits formed on the prism (1010) faces of ice needles were 
rectangular in cross-section and mostly flat at the bottom as shown in 
fig. 9. Star-shaped etch channels were seen only very rarely in crystals 
grown from the vapour in virtually dust-free conditions; one such is 
shown in fig. 10. The channels run parallel to the (1120) axes of the crystal 
and are marked all along their length by very small pits less than 1 in 
diameter. 

In an attempt to correlate the positions of etch pits arising on opposite 
sides of a plate crystal, about 100, thick, the crystal was suspended from 
a glass fibre and coated on both sides with Formvar solution. The plastic 
film was then removed by folding transparent adhesive tape round the 
crystal, peeling it off, and straightening the tape to obtain replicas of 
both faces side by side. Although both basal faces showed approximately 
the same density of etch pits, it was not possible to match their positions. 
This suggests that, if the pits were formed by the etching of dislocations, 
the latter followed rather devious paths in traversing the 100» thickness 
of crystal. 


§5. THe MoveMENT OF DISLOCATIONS UNDER MECHANICAL STRESS 


Replicas of the surfaces of thin sheets ofice, 1-2mm in thickness, produced 
by the freezing of water in a refrigerator, revealed large concentrations 
(~ 105/cm?2) of etch channels oriented along the (1120) directions as shown 
in fig. 11. The specimens were subjected, at —20°c, to quite strong 
vibrations from the compressor of the refrigerator ; when similar specimens 
were mounted on foam rubber, so that the vibrations were greatly damped, 
the number of etch channels was greatly reduced. We believe that the 
channels, which were tapered along their length, resulted from the move- 
ment of edge dislocations under mechanical stresses induced by the 
vibrations. During a period of 20min the etch pit advanced about 15 
during which time the specimen was subjected to about 30 000 vibrations. 
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The average movement per vibration was therefore about 5A, or about 
one inter-molecular spacing along the <1120) direction! 

Similar sheets of ice were subjected, again at — 20°c, to local stresses 
by dropping lead shot onto the surface from heights of a few inches. The 
replicas now revealed rosettes of the type shown in fig. 12 and noted earlier 
by Gilman and Johnston (1957) on the surfaces of lithium fluoride crystals. 
Arrays of etch channels, presumably associated with similar arrays of 
dislocations, radiated from large etch pits along the (1120) directions which 
again appear to be preferred glide directions. 


§ 6. Discussion 


The fact that the surface densities of the smaller etch pits on ice lie 
in the range 10° to 108/cm? and are comparable to the observed densities 
of dislocations in other materials, suggests that they probably arose 
by the etching of emergent dislocations. A close association between 
dislocations and etch pits is also suggested by the latter tending to group 
preferentially along grain boundaries and slip lines. 

The inverted pyramidal shape of the pits on the basal plane may be 
interpreted in terms of a dislocation being etched at similar rates in direc- 
tions normal and parallel to the surface, the dislocation following a line 
drawn from the apex of the pyramid to the centre of the base. Anisotropic 
lateral etching on the prism faces is revealed by the rectangular shape of 
the pits. 

In the present work it was usually not possible to determine whether the 
pits were generated about edge or screw dislocations but etching patterns 
sometimes appeared which indicated the presence of helical dislocations. 
These were first discovered by Bontinck and Amelinckx (1957), their 
formation being explained by Amelinckx, Bontinck, Dekeyser and Seitz 
(1957) who suggested that the condensation of vacancies on to a dislocation 
possessing a strong screw component could cause it to climb into a helix. 
Etching of such a dislocation, which may be regarded as an unwinding of 
the screw component, would then produce a pit in the form of a helical 
cone. Amelinckx, Bontinck and Dekeyser (1957) further demonstrate that 
bunching of the unwinding monolayers, caused by the rotation of the centre 
of etching along the helix, could form visible spiral steps whose spacing 
would be equal to the pitch of the helix and, in consequence, bear no direct 
relation to the Burgers vector of the underlying dislocation. This mechan- 
ism may well explain the appearance of the spiral steps in fig. 3, the height 
of which were of order |, and therefore comparable to the pitch of helical 
dislocations observed directly in other crystals. 

It is also known that two parallel helices of opposite sign may combine 
to produce a row of closed loops and should this emerge from the etched 
surface at a slight angle to the vertical, it could give rise to etch pits con- 
sisting of a series of closed terraces as shown in fig. 6. 

The double rows of symmetrical etch pits seen in fig. 10 may have 
originated on helical dislocations with axes parallel to the surface. In this 
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and similar cases they may have been initiated by a dust particle in much 
the same way as Jones and Mitchell (1958) were able to generate helical 
dislocations and rows of prismatic loops by pressing small glass spheres 
into crystals of silver chloride. 

Our attempts to reveal the three-dimensional structure of dislocations in 
ice, by adding to the water small quantities of silver nitrate in the hope 
that the silver ions might migrate preferentially to the dislocations line 
during the freezing process, have been unsuccessful. 

The pyramidal hillocks, which appeared in concentrations of between 
10? and 10°/cm? on the basal surfaces of crystals grown from the vapour 
in very clean conditions, were not observed on crystals grown from the 
melt. Although they may have been formed by etching during the 
preparation of the first replica, the fact that they were eroded and finally 
removed by subsequent etching leads us to think that they may have 
been formed by growth rather than by dissolution at preferred sites on 
the crystal surface. This point is not easy to check but is being investi- 
gated further. 

The fact that dislocations spreading out from localized regions of high 
stress and etch pits moving under applied stresses show a very strong 
tendency to follow the (1120) directions suggests that these are preferred 
glide directions. 
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ABSTRACT 


Measurements have been made of the variations with temperature of the 
intensities of x-ray diffraction lines from silver and iron. For silver the 
measurements were in the range 286—-1100°K, and the Debye—Waller 
temperature factor corresponded to a characteristic temperature of 197°K at 
286°xK. For iron the measurements covered the range 286-1190°K, and 
corresponded to © =389°x, although there was a large scatter in the results 
due to crystal changes brought about by annealing. The intensity measure- 
ments showed no discontinuity at the Curie temperature, nor was there any 
anomaly in the spacing of the (110) planes. 


§ 1. INTRODUCTION 


THE variations with temperature of the integrated intensity and peak 
height of x-ray diffraction lines 331, 551/333 and 422 for a silver powder 
specimen and 220 for an iron specimen, were measured at high tempera- 
tures. For silver an explicit volume dependence of the characteristic 
temperature was adequate to give the observed temperature dependence 
in terms of the Debye-Waller theory, the value © = 197°K being obtained 
at room temperature. This result does not agree with the work of 
Andriessen whose results give © = 239°K when recalculated to allow for 
lattice expansion, nor with the work of Boskovits et al. (1958) who 
obtained a large deviation from the theoretical intensity variation above 
700°k. The present work does agree approximately with that of 
Spreadborough and Christian (1959 b). 

The results for iron showed a much larger scatter than silver, associated 
with crystal changes brought about by prolonged annealing. No anomaly 
was observed at the Curie temperature, either in the intensity variations 
or in the spacing of the 110 planes. 


§ 2. THEORY 


For the Debye-Waller theory, the temperature dependence of the 
intensity of an x-ray reflection is given by the factor exp (— 2) (following 


James 1954, 5.67), where 
6h?2T sin? 6 x 
by mk Oxy? ee . {te + i} . . . . . . (1) 
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The characteristic temperature @,, has been regarded as a constant in 
the above expression: amongst other things this is equivalent to assuming 
that the volume of the crystal remains constant as the temperature 
changes. Paskin (1957) has given a thermo-dynamical derivation of the 
variation of @,, (regarded as a temperature dependent parameter @(T)) 
with volume, the relation for an isotropic crystal which obeys Gruneisen’s 
expansion equation being written 
O7) = (22)" : in) 
(75) arp ; 
where a, is the lattice parameter at temperature 7°K. The Gruneisen 
constant y (following Hume-Rothery 1945) is calculated from the equation 
VB 
Lem ey? 
V being the atomic volume, CO, the atomic specific heat at constant value, 
B the volume coefficient of expansion and k the compressibility. The 
values used for the constant y are 2-40 for silver and 1-6 for iron. 
Substituting for ©,,? in eqn. (1) we have 
6/27 sin? @ : 

letra es et rena 
where 7” =7'(a)/a7,)*" the ‘adjusted’ temperature allowing for the effect 
of lattice expansion, and 

Dix} = (p(w) + w/4} 


is approximately unity and is tabulated by James (1954). 

Neglecting the small effect due to change of Bragg angle @ with tem- 
perature, the intensity p, of a given reflection at temperature 7' can be 
compared with the intensity at a temperature 7'o, the relation being 

pr\__—«‘12h? sin? 0 Digy .s : 
in (2) boat fe | sey Te ik Nei 2 
Thus a linear relation is expected between the parameter 


A® PS 
Y= —— — 
sin” 0 h (S 


and the variable 7" = (O{x}/®{x)}7". A very approximate value of @,, is 
satisfactory for calculating values of ®{z}, since this function is nearly 
unity for all values of x corresponding to temperatures above ®,,. Any 
departure from a linear relation can be attributed to an additional 
dependence of ©,, on temperature, above that related to the coefficient of 


expansion. From the slope of a linear relation a value of @, can be 
obtained. 


§ 3. EXPERIMENTAL MEASUREMENTS AND RESULTS 
3.1. Silver 


Filings were prepared from pure (99-99%) silver, carefully cleaned and 
that portion which passed through a 350 mesh sieve prepared into a flat 
diffractometer specimen, held in a molybdenum trough. 
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This was mounted in a furnace built onto a Norelco diffractometer, 
described by Spreadborough and Christian (1959 a), and the specimen was 
annealed at 1020°K for ten hours in vacuo. It was found that subsequent 
annealing treatments at temperatures up to 870°K did not alter the room 
temperature line profiles of the 331, 422 or 511/333 reflections obtained 
using CuKa or MoKa radiation. 

The integrated intensity for a given reflection obtained for a specimen 
temperature 7’) could then be compared with that obtained for a higher 
temperature 7’, using a planimeter to measure the total area under the 
4%, profile, above the background level in the immediate neighbourhood 
of the reflection. From these measurements values of Y corresponding to 
different temperatures could be obtained. 

The background intensity in the neighbourhood of the line was observed 
to increase with increasing temperature of the specimen, but the variation 
of intensity with angle in the ‘wings’ of the Bragg peak was small. 
However, no attempt was made to obtain an absolute value for the Bragg 
scattering, but the variation with temperature of the intensity, measured 
over a small fixed angular range, was found. 

The error due to including temperature diffuse scattering in the measure- 
ment of the integrated intensity has been discussed by Chipman and 
Paskin (1959), who showed that the error introduced by including all 
the wings of the lines (a procedure described by Renninger (1952)) would 
lead to a larger value of the characteristic temperature than the true one. 
However, if the temperature diffuse scattering over the small angular 
range measured is calculated, then the error in characteristic temperature 
caused by including all the temperature diffuse scattering in the measure- 
ments is less than 1°% for any of the Bragg reflections used. For the 
experiments where the peak heights were compared, the error could be 
larger. The correct adjustment of the apparatus and the validity of 
background level was checked by observing the 2:1 intensity ratio of the 
a2, doublet lines. 

Values of Y could also be calculated from variation in peak heights with 
temperature. The comparison between measurements based on inte- 
grated intensity values and on peak heights was made because : 

(a) The theoretical calculations refer to the integrated values. 

(b) Peak height measurements can be made more accurately than a 
complete profile, because small fluctuations in temperature (~ + 2°C) 
during the profile measurements cause change of Bragg angle which 
result in serious errors. Provided surface effects in the specimen did 
not produce asymmetry in an x-ray line, then instrumental effects are 
unlikely to introduce errors when comparing peak heights instead of 
integrated intensities. However, it was found that changes in crystalline 
arrangement introduced by annealing treatments produced greater changes 
in peak heights than in integrated intensities, and great care was necessary 
in making use of the peak height data. After repeated annealing at 
1220°K the measurements were successfully extended up to 1100°K. 
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The graph of Y against the adjusted temperature T” is shown in fig. 1, 
and is seen to be linear up to a specimen temperature of 1100 °K, the highest 
temperature reached. The slope of the linear relation corresponds to a 
characteristic temperature of 197°K + 6 at 286°K, whilst the usually accepted 
specific heat value is 210°K (James 1954). Other x-ray values for the 
characteristic temperature of silver, taken from the literature, but re- 
calculated using the method of analysis used in this paper, are: 

Andriessen (1935) Oy = 239°K + 5, 

Boskovits et al. (1958) Oy, = 188°K + 12, 
using only data for temperatures above 290°K, and 

Spreadborough and Christian (1959b) ©=191°K £8. 


Fig. 1 
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Of these results that of Andriessen gives a nearly linear relation up 
to 880°K, but the characteristic temperature obtained is much higher than 
expected. The work of Boskovits et al. gives an approximately linear 
relation from 290°K to 680°K, the slope of the relation agreeing with the 
present work to within the experimental error. Above 680°K, however, 
measurements at 780°K showed a more rapid diminution of x-ray intensity 
than expected, but such a derivation from linearity is not supported by the 
present work. 

The method of correcting for the background used in the present work 
would tend to give too high a value for @,,. The value of the characteristic 
temperature used in X-ray measurements is expected to be larger than that 
used in the Debye specific heat theory, owing to different methods of 
averaging the velocities of longitudinal and transverse elastic waves in the 
crystal: in the case of silver we expect the X-ray characteristic temperature 
to be 4% greater than that used in specific heat measurements, from the 
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work of Zener and Bilinsky (1936). The observed characteristic tem- 
perature pertaining to the x-ray temperature factor is, however, lower 
than the specific heat value. 


3.2. Iron 

Two iron specimens were prepared and used: one from Johnson Matthey 
iron sponge having an analysis 0-0053% Ni, 0-0003°%% Mn, 0-063% insoluble, 
and the other by filing B.I.8.R.A. H-iron under cold paraffin. In both 
cases the powders were carefully cleaned and then passed through a 350 
mesh sieve before preparing into a specimen by pressing into a pure 
iron trough. The second specimen was vacuum-annealed for four hours 
at 1400°c before placing in the x-ray camera. 

Using filtered iron radiation, the 220 reflection at approximately 145° 20, 
was measured. The integrated intensities at different temperatures were 
compared, the measurements being taken from the trace on an automatic 
recorder chart, corrected for counter dead-time losses. The results showed 
a larger scatter than in the case of silver, and a large number of chart 
recordings were used: 


6» 328°K 


4 Fe specimens. 
Fe radiation. 220 reflection. 


A graph of parameter Y against ‘adjusted’ temperature 7” was plotted, 
as in the case of silver, the adjustments to the temperature scale being 
made using the lattice parameter values of Basinski e¢ al. (1955) y= 1-60 
and ©=430°K. The value of the characteristic temperature obtained 
from this graph would depend only to the second order upon the © value 
used in making the temperature adjustment. 

The scatter of the results, as presented on this graph is seen to be much 
larger than in the case of silver, but does show that there is no marked 
anomaly in the temperature factor at the Curie temperature. This large 
scatter could have been due to statistical changes in the distribution of 
reflecting particles, but there was also a time-dependent, partly reversible 
effect, superimposed on the random variations. This is shown by the 
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measurements represented by crosses on the graph, which were made 
immediately after heating to a higher temperature, whilst the other points 
correspond to measurements made after an annealing period of one hour 
or longer. This effect may be due to extinction variations brought about 
by polygonization or other changes in the iron crystals. 

Fitting a straight line passing through 7',=286°K to these results, and 
using the same analysis as for silver gives a ©,, value for the points repre 
sented by crosses of 389°K + 17, and for the other points @,,=328°K + 20. 
The largest characteristic temperature which can be fitted to any 
measurement is 403°K, whilst the accepted specific heat value for «-iron is 
430°K. It would seem likely that the discrepancy may be associated with 
changes in extinction brought about by annealing. As in the case of 
silver, inclusion of temperature diffuse scattering in the measurements 
would lead to too high a value of ©,, ; although the error due to including all 
the diffuse scattering, calculated following Chipman and Paskin (but for a 
b.c.c. lattice) over the range measured of 1-6° 26, would lead to an error in 
©,, of less than 0:3%. 

The variation in Bragg angle of the a, and a, peaks with temperature 
gave a measure of the thermal expansion of the 110 direction in the iron 
crystal, with a sensitivity of +0-0002kx in d value. No anomaly in 
thermal expansion was observed at the Curie temperature, in agreement 
with the work of Basinki et al. (1955), who used an extrapolation technique 
such that their lattice parameter results were effectively based upon a study 
of the 310 direction. 
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ABSTRACT 


In a crystalline element the scattering of long wavelength neutrons by a 
point defect and an atom are known to be identical. However, when lattice 
distortion around the point defect is also considered, this is no longer true; in 
particular the scattering of a point defect is now anisotropic and varies with 
wavelength. An expression for the cross section of a defect in a relaxed 
lattice has been derived and applied to f.c.c. and b.c.c. lattices, and in greater 
detail to copper and molybdenum where the magnitude of the relaxation 
around defects is known. It has been found that the scattering cross section 
is appreciably altered by the relaxation of the surrounding atoms, and in 
particular from those that are nearest neighbours to the defect.. A neglect of 
this relaxation in the interpretation of a neutron scattering experiment could 
lead to an error of an order of magnitude in an estimate of the number of 
defects. Also a study of the variation in cross section with wavelength or 
scattering angle, which is analogous to an x-ray small angle scattering experi- 
ment, is incapable of determining uniquely the size of defect agglomerates in 
solids. However, under favourable circumstances, such experiments could 
become a useful technique for estimating the lattice distortion around defects. 


§ 1. INTRODUCTION 


Many experiments on solids irradiated by energetic particles have aimed 
at estimating the numbers and states of agglomeration of the resulting 
point defects. However, a unique description from experiments observing 
changes in physical properties has proved very difficult; probably the 
most unequivocal have been observations under the electron microscope, 
but unfortunately they are limited to agglomerates larger than about 
154. 

The scattering of long wavelength neutrons, for which Bragg reflections 
are absent, is a particularly attractive tool for studying smaller 
agglomerates and elementary point defects in crystalline solids, because 
the scattering cross section per point defect is the experimentally known 
bound atom coherent scattering cross section. Experimentally the main 
difficulty in the past has been the low flux of long wavelength neutrons 
available from a reactor; this technique is likely to be exploited more, 
however, as cooled moderators (Butterworth e¢ al. 1957) inside high flux 
reactors become available as neutron sources. 

Initial experiments by Antal e al. (1955) were on neutron irradiated 
graphite. Their estimate of the defect concentration was in good 
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agreement with theoretical predictions. They pointed out, however, 
that it may be in error for three reasons: (a) irradiation may alter the 
inelastic scattering cross section of the solid, (b) defects may be aggregated 
in pairs or larger defects, (c) there may be relaxation of atoms in the 
neighbourhood of the defect. Changes in the inelastic scattering cross 
section have been considered negligible because of the small change in 
Debye temperature after irradiation. Antal and Goland (1958) have 
discussed the effect of defect pairs on the wavelength dependence of the 
transmission cross section; they were able to determine the separation 
of defect pairs in neutron irradiated a-Al,O3, although the complexity 
of the material was an additional difficulty in the interpretation. 

We propose to discuss here the third possible source of error, namely 
the effect of lattice relaxation about the defect on the scattering cross 
section. Unfortunately the magnitude of this relaxation is uncertain 
in those materials that have been studied experimentally. However, 
it is, known approximately in certain metals, and for two of these, copper 
and molybdenum, the cross sections are calculated in this paper to 
illustrate the general principles. 


§ 2. THEORY 


Let us consider a perfect lattice of N atoms per unit volume to which 
v atoms have been introduced randomly on interstitial sites, causing 
relaxation of » atoms around each interstitial. We note here that the 
ensuing treatment is identical for v vacancies. 

The scattering amplitude per unit solid angle for unit incident flux 


Ni v 
=o > exp(—ix.e,)+ > «SER loa iba) 


B=1 q=1r= 
- » Sexpl-% E, +901} . eee) 
= b[4+B-C] ee Peds 


where 6 is the bound atom scattering amplitude per unit solid angle for 
unit incident flux, x is the scattering vector, magnitude (47sin6)/A 
(20=scattering angle, A=neutron wavelength), pe, and &, are vector 
distances from the origin of atoms in lattice and interstitial sites 
respectively, tp, is the vector distance of a relaxed atom from its neigh- 
bouring interstitial and «, is the corresponding distance of a lattice site 
vacated due to relaxation. The sum over p is for the perfect, unrelaxed 
lattice ; that over r includes a term with w =0, which corresponds to the 
contribution made by the interstitials to the scattering amplitude. 
The scattering intensity is 


6°|AA* + A(B*—C*)+ A*(B—O) + BB* +CC*— BC*—CB*]. (2) 
Beyond the Bragg cut off wavelength 
AA*~1, 
A(B*—C*) + A*(B—C) x v2, 
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and turn out to be negligible compared with the other terms, which may 
be expressed thus: 


b > exp[—ix.€,-8,) S exp[—ie.(h,-o,)] 


aq =1 rr’ =0 
+ ¥ expl-ie(p.-@l- 5 exp[—ix. (b,-@,)] 
8, 8’= r=0,s=1 
— _S expl-%.(@-bal} 5 ee) 


The sum outside the brackets is of v? vectors of random phase, apart 
from the restriction that §,, §,, must lie on interstitial sites; it may be 
shown to equal [v(N —v)]/N. Hence in most practical cases, when vy is 
less than a few per cent of V, the sum equals v to a good approximation, 
which is the value for a perfectly random array of v interstitials. 

The four series in brackets, which consist of three sums of the general 
type cos [x.(x—y)], are correct for a single crystal, when the orientation 
of the relaxed zone about each interstitial is identical. For polycrystalline 
material we must average over all orientations of (x—y), resulting in a 
series of terms of the. type 

sin |||. |x —y|] 
be] ey] 

Hence the scattering intensity per unit incident flux, per unit solid 
angle, which is the differential scattering cross section is 


. (4rsin 8 . (4rsin @ 
[ sin( = lb, i v,| ) fa <in( = x lep, ze el) 


o, = bv Sy ae oe > > 
Ee a kb, b,-| Laie Ss lp;— Pol 
. f478in8 
/ <in( He \v.— eel) 
LM Dene snmp 5 sao ee (4) 
47 sin 0 
tier trsinty, ol 


The moduli represent distances within the relaxed region of the crystal 
surrounding one interstitial, and the three sums correspond to all 
distances between atoms, lattice sites, and an atom and lattice site 
respectively. 

Integrating each term of (4) over 47 of solid angle we obtain the total 
eross section, which is obtained from a transmission experiment : 


/2 A 2 
sin F lee ) ia sint(Z \e.— el) 
7,7 =0 & lve ) , + les— ¥el 
simn(F Ive) 
Sp Sf al TAS 
an (= Il) 


—2 


The sums are identical to those in (4). 


P.M. sm 
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These summations are similar to those in the diffraction theory of 
polyatomic molecules (e.g. James 1948), where terms corresponding to 
all interatomic distances in the molecule are included; in comparison 
the sums here include vacant lattice sites, and some terms are negative. 

Finally we observe that in a lattice containing two sorts of defects, 
both with relaxed neighbouring atoms, for example interstitials and 
vacancies, the total scattering is simply the sum of two series such as 
(4) or (5) to a good approximation, provided that the defect concentration 
is not greater than a few per cent. 


§ 3. APPLICATION TO FacE-CENTRED AND BobY-CENTRED CuBIC MuTats 

We present here values of the total and differential scattering cross 
sections for a point defect in a f.c.c. and a b.c.c. metal considering 
relaxation of nearest neighbours only, the variable parameters being the 
relaxation «, the wavelength and the scattering angle. 

Figure 1 depicts the total cross section of an interstitial in copper as a 
function of neutron wavelength. Three values of «, the distance that the 
six nearest neighbours have moved outwards from their lattice sites, 
are shown; in all cases « is expressed in units of half the lattice spacing, a. 
For copper Huntington (1953) and Tewordt (1958) find values of « ranging 
from 0-210 to 0-257. Cross sections are expressed as fractions of the very 
long wavelength value 47b?; this is also the value o, that would be 
observed if the lattice was not relaxed. Values of a/A on the abscissa 
enable these results to be applied directly to other f.c.c. lattices. 

Figure 2 shows the variation in total cross section with « for three 
values of a/A. 

The solid lines in fig. 3 depict the differential scattering cross section 
of an interstitial in copper as a function of angle. They are for a neutron 
wavelength of 7-6A, which corresponds to a/A=0-475. 

Figures 4-6 are the corresponding curves for a vacancy in copper, 
where the twelve nearest neighbours are relaxed inwards. Tewordt 
(1958) has found e=0-0194/2 and 0-0234/2, using two different potentials. 

Figures 7-9 are the corresponding curves for a vacancy in a b.c.c. 
material with the eight nearest neighbours relaxed and A and @ values 
that are applicable to molybdenum. The values of « may be compared 
with that of 0-05054/3 calculated by Girifalco and Streetman (1958). 


§ 4. Discussion 


The results of figs. 1-9 show how marked is the influence of reasonably 
expected values of relaxation on the scattering cross section. A neglect 
of it in metals could lead to an experimental estimate in the numbers of 
defects that was in error by a large factor. 

It is pertinent to ask how much in error these results are compared with 
physical reality, when other neighbouring atoms are relaxed as well. 
Let us consider the 6(1,0,0) and 24(2, 1,0) neighbours about an inter- 
stitial in copper relaxed to new positions (1+a, 0, 0), (2+, 1+7, 0), 
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where coordinates, centred about the interstitial, are in units of a/2. 
Tewordt (1958) has shown that these are the two kinds of atoms most 
relaxed. Transmission and differential cross sections have been calculated 
for certain values of \ and 0, for «=0-25, w=0-050, 7=0-055, which are 


Fig. 1 


ovis) 
ot 
4b? 
eS) 
@nS 


5 
\ for cu,& 


So} 
0.4 0.5 0.6 0.7 O23 ao) 


The total cross section o, of an interstitial in a f.c.c. lattice plotted as a function 
of wavelength. , is divided by 47b?, which is the total cross section of 
an interstitial in an unrelaxed lattice. The curves are results of cal- 
culations when only nearest neighbours are allowed to relax, The 
points make allowance for the relaxation of other neighbouring atoms as 
well, when «=0:25. Values of these relaxations which are valid for 
copper were used. 
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values lying between those calculated by Tewordt using two different 
potentials. The results are presented as points on figs. | and 3, and 
may be compared with the curves for «= 0-25, when relaxation of (2, 1, 0) 


atoms is neglected. 
Fig. 2 
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ioe ee has also calculated the relaxation of next-nearest neighbours 
ag a vacancy In copper. The nearest and next nearest neighbours 
coordinates (1,1,0) and (2,0,0) are relaxed to (1—B8, 1-8, 0) an 
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(2—€, 0, 0) respectively, where 8=0-020 and £=0-002. The results are 
shown on figs. 4 and 6 and are to be compared with €=0-02,4/2, when 
relaxation of (2,0,0) atoms is neglected. 

For a vacancy in molybdenum the coordinates of the relaxed nearest 
and next-nearest neighbours are (1 —y,l—y,l1—y) and (2—y,0,0) 


Fig. 3 
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The differential cross section o, of an interstitial in copper, plotted as a function 
of angle. o, is divided by b?, which is the differential cross section of an 
interstitial in the unrelaxed lattice. The neutron wavelength is 7-6 A. 
Solid lines are results of calculations when only nearest neighbours are 
allowed to relax. The points make allowance for the relaxation of other 
neighbouring atoms as well, with «=0-25. The broken line is for an 
interstitial in a copper single crystal for neutrons incident in the [100] 
direction and scattered in the (001) plane. Relaxation of nearest 
neighbours only is considered. 
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respectively, where y= 0-0505 and y= 0-0224, using the results of Girifalco 
and Streetman (1958). The results may be compared with the nearest 
neighbour approximation for «=0-05,/3. 

These results indicate that the cross section is only slightly in error by 
neglecting the relaxation of next nearest neighbours. Only certain parts 
of the curves for an interstitial in copper, when the relaxation is very 
large, are measurably affected. 


Fig. 4 
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The total cross section o, of a vacancy in a f.c.c. lattice plotted as a function of 
wavelength. The solid line is the result of calculations when only 
nearest neighbours are allowed to relax. The points make allowance 
for the relaxation of next-nearest neighbours also, using values which 
are valid for copper. For comparison the broken line shows the total 
cross section of a vacancy pair in a f.c.c. lattice, assuming there is no 
relaxation of the surrounding lattice. 


We do not expect the relaxation of atoms a large distance from the point 
defect to influence these scattering cross sections. This is because in 
the elastic continuum approximation the density of an infinite medium 
at any point remains unchanged when a point defect is introduced into 
it, so that the scattering intensity is zero, as in the perfect lattice. 

All the results presented so far have been for randomly oriented 
polycrystalline solids. To evaluate the differential scattering cross 
section for a single crystal we must sum terms of the general type 
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The total cross section of a vacancy in a f.c.c. lattice plotted as a function of the 
relaxation. The calculations make allowance for the relaxation of 
nearest neighbours only. For copper «~ 0-028. 
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The differential cross section o, of a vacancy in copper, plotted as a function of 
angle. The neutron wavelength is 7-6 4. Solid lines are results of cal- 
culations when only nearest neighbours are allowed to relax. The points 
make allowance for the relaxation of the next-nearest neighbours as well, 
with «=0-024/2. For comparison the broken line shows the differential 
cross section of a vacancy pair in copper, assuming there is no relaxation 
of the surrounding lattice. 
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a 


cos [«.(x—y)]. The dashed curve in fig. 3 illustrates the differential cross. 
section of an interstitial in a copper single crystal for -=0-25. The beam 
is incident along the [100] axis and is scattered in the (001) plane. It is 
seen that for small phase differences, when 6 is small, the scattering 1s 
the same as polycrystalline material; only at larger values of 6 do the 
curves diverge. 


Fig. 7 
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The total cross section o, of a vacancy ina b.c.c. lattice, plotted as a function of 
wavelength. The curves are results of calculations when only nearest 
neighbours are allowed to relax. The points make allowance for the 
relaxation of other neighbouring atoms, when «=0-054/3. Values of 
these relaxations which are valid for molybdenum were used. 


It was pointed out above that the cross section of agglomerated defects 
will depend on wavelength and angle. The scattering from a defect pair 
(neglecting relaxation) in copper is shown as a function of \ and @ in figs. 
4 and 6 respectively, for comparison with relaxation effects. Apart from 
low scattering angles the curves are quite different, so that in principle 
the two effects are experimentally separable (but see next paragraph). 

The results of this paper show that lattice relaxation and the state of 
agglomeration of defects influence the scattering cross section to a 
comparable extent. This is a serious disadvantage when neutron 


oO 
0-02 0-04 0:06 0.08 0.10 € 


The total cross section of a vacancy in a b.c.c. lattice plotted as a function of the 
relaxation. The calculations make allowance for the relaxation of 
nearest neighbours only. For molybdenum «~ 0-087. 


Fig. 9 
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The differential cross section o, of a vacancy in molybdenum plotted as a 
function of angle. The neutron wavelength is 7:64. Solid lines are 
results of calculations when only nearest neighbours are allowed to relax. 
The points make allowance for the relaxation of the next-nearest neigh- 
bours as well, with «=0-054/3. 
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scattering experiments are used to study the nature of defects produced 
after irradiation damage because these two effects cannot be distinguished, 
since only for certain elementary defects in a limited number of materials 
is the relaxation known. To illustrate this, fig. 4 depicts the scattering 
of a vacancy pair, without relaxation, the value of which is unknown. It 
may be that with relaxation included the curve would become a straight 
line independent of wavelength, like a defect with no relaxation, or 
alternatively have the shape like the other curve in fig. 4 which is for a 
single vacancy with relaxation, so that results would be ambiguous. 
Perhaps in certain particular cases, for example in materials possessing 
a more open structure, when relaxation should be less pronounced, or 
with single crystals, where the cross section fluctuates more rapidly with 
wavelength and scattering angle there may be less ambiguity. 

However, if the size of defect agglomerates is known in a particular 
solid, measurements of the neutron scattering cross section will yield 
information about the relaxation of the surrounding lattice. One such 
example is a very dilute solid solution, when the solute atoms are effectively 
point singularities randomly dispersed in the solvent lattice ; in particular, 
if this lattice possesses cubic symmetry, the relaxation of nearest neigh- 
bours around each solute atom may be determined directly from such an 
experiment. 


REFERENCES 


AnTaL, J. J., and Gonanp, A. N., 1958, Phys. Rev., 112, 103. 

AntaL, J. J., Weiss, R. J., and Dimrnzs, G. J., 1955, Phys. Rev., $9, 1081. 

Burrerworts, I., Eceistarr, P. A., Lonpon, H., and Wess, F. J., 1957, Phil. 
Mag., 2, 917. 

GrriFaLco, L. A., and SrrEErMAN, J. R., 1958, J. Phys. Chem. Solids, 4, 182. 

Huntrineron, H. B., 1953, Phys. Rev., 91, 1092. 

ieee Es The Ortical Principles of the Diffraction of X-rays (London: 

ell), p. 482. 
Teworpt, L., 1958, Phys. Rev., 109, 61. 


eles 


The Nature of Radiation-induced Point Defect Clusterst 


By R. S. Barnzs and D. J. Mazzy 
Atomic Energy Research Establishment, Harwell, Didcot, Berks 


[Received June 21, 1960} 


ABSTRAOT 


The electron microscope has been used to examine point defect clusters 
produced in copper and aluminium foils bombarded with 1-4 x 1017 alpha- 
particles em-?. The foils were bombarded as a stack and, in copper, both 
those through which the alpha-particles had passed and that in which they 
came to rest appeared similar, containing dislocation loops (r~200 A) and 
a background of more numerous and smaller dots (r~20 4). The behaviour 
at grain boundaries suggested that these two types of cluster were due to 
different point defects. Whereas the dislocation loops normally annealed 
out at about 350°c, in the foil in which the alpha-particles came to rest 
(containing ~10! atoms of helium cm~—?) they grew to form a dislocation 
tangle, and eventually small helium bubbles (r~40 A) appeared in the same 
number as the original small dots. These results indicate that the dislocation 
loops are due to the clustering of interstitial atoms and the small dots are 
clusters formed from vacancies. It is also inferred that the helium 
atoms are in interstitial positions during the bombardment, and form 
bubbles by nucleating upon the clusters. 


§ 1. IyTRODUCTION 


THE much smaller energy needed to form a vacancy ensures that the point 
defect clusters observed in quenched metals are due to precipitation of the 
excess concentration of vacancies and not interstitial atoms. However, 
during atomic bombardment vacancies and interstitial atoms must be 
produced in nearly equal numbers, because they are produced by a collision 
process, and the effects observed may be due to either type of defect. 
Bombardment with alpha-particles introduces both point defect clusters 
and helium atoms, the latter on subsequent heating capture large numbers 
of vacancies (Barnes et al. 1958) and enable the nature of the point defect 
clusters to be determined. 

Thin films of copper extracted from a solid block which had been 
bombarded with alpha-particles have been examined and helium bubbles 
observed within them (Barnes 1960). However, the extraction of thin 
films from a solid block is tedious and wasteful of material. The results 
described here were obtained by placing stacks of thin sheets of metal in the 
external beam of a cyclotron accelerating alpha-particles. These were 
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then thinned by normal electrolytic methods until they could be viewed in 
transmission in the electron microscope. Sheets from various depths 
beneath the bombarded surface were studied and the effects produced by 
the passage of the alpha-particles (in the sheets near the surface) compared. 
with those produced where the helium atoms were deposited (in the sheet 
corresponding to the end of the range of the alpha-particles). 

Johnson Matthey spectroscopically standardized copper or aluminium 
was used and the sheets were initially either 0-002 in. or 0-001 in. thick and 
were vacuum annealed for 2 hours at 650°c before bombardment; during 
bombardment their temperature did not exceed 250°c. As the range of the 
38 Mev alpha-particles from the cylotron was 0-009in. in copper, about 
five 0-002 in. or nine 0-001 in. sheets were affected during each bombardment. 

However, normally only one of the sheets contained the helium. This 
was located by observing the helium bubbles after heating a little of a 
sheet to about 800°c for 2 hours. Figure 1+ shows the cross section of a 
0-001in. foil (after electrolytic polishing) containing unresolved helium 
bubbles near the grain boundaries (Barnes et al. 1958). Suitable masking 
enabled other pieces from this same sheet to be thinned (after the appropriate 
annealing treatment) leaving the centre of the helium layer within it. 


§ 2. RESULTS 

The photographs in figs. 2 and 3 were taken from foils, containing no 
helium, through which 1-4 x 10!” alpha-particles em? (32 4 amp hours) had 
passed. The appearance of the electron micrographs taken from various 
distances along the alpha-particle tracks showed little variation, each 
containing a number of relatively large dislocation loops (7 ~ 200 A) and a 
background of much smaller dark dots (7 ~ 20 A) (fig. 2). Figures 3 (a), (6) 
and (c) are successive photographs of a similar foil taken while it was at 
about 350°c in the ‘hot-stage’ of the electron microscope. At this tem- 
perature the loops disappeared in less than 30 min. Sheets were also 
annealed in a vacuum furnace at various temperatures. After a 2-hour 
anneal at 300°c there was little change in the appearance of the films, after 
the same time at 325°c slight annealing had occurred and after 4 hours 
at 350°c complete annealing of the loops had occurred. 

The thin films in which the helium atoms had come to rest contained both 
large loops and small dots, which were not significantly different in size 
from the above (table), but the dots appeared to be localized around the 
loops (fig. 4). The annealing behaviour was however markedly different 
and the series of photomicrographs of foils annealed in a vacuum furnace 
shows this (fig. 5). Instead of the dislocation loops disappearing after 
2 hours at 350°c (fig. 5a) they enlarged and after 2 hours at 450°c (fig. 5b) 
had become a tangle of dislocations. After this anneal there were small 
bubbles near the grain boundaries (fig. 5b). After 2 hours at 550°c small 
bubbles pervaded the grains and on further heating for 2 hours at 750°c they 
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became larger (fig.5c andd). At this temperature larger bubbles occurred 
along the grain boundaries (fig. 5e) consistent with previous observation 
(Barnes 1960). 

Aluminium foils through which the alpha-particles had passed contained 
no observable damage effects, but the foil with helium in, contained a 
dislocation tangle, and upon heating small helium bubbles appeared. 


Number of atom spaces 
Radius Number Ming tog 
elects a) Loops Spheres 
(477//3)(r/a)?n| (1677/3) (7/a)8n 
Point 
Figure 2 (a) | defects | 2:-2x10-6 | 4:2 10° Pet 10" — 
only 
Point 
Figure 2 (6) | defects | 2-2x10-* | 6-4x10® | 1:7x 104 = 
only | 
He and 
Figure 4 point LG 5e eee teed xe Ae e201 04 — 
defects 
Point 
Higure 2(a): detects’ | 2-5 10>?) 79-0 «101 23-610" 5-510 
only 
He and 
Figure 4 point | I-9x 10" | V-2x104 1 2-4 x10! a0 10 
defects 
Figure 5 (d)|, #® | 4-0x10-7 | 14x10" se 3-2x 1015 
igure © (¢) | bubbles 


The size and number per unit area of the dislocation loops (the first three rows), 
the dots (the next two rows) and the bubbles (last row) observed in the thin 
films. The foil thicknesses were approximately 1000 A. 


§ 3. Discussion 


Dislocation loops form in many metals during quenching (Hirsch et al. 
1958) their number and size depending upon the quenching conditions. 
These loops which are produced by the precipitation of the highly super- 
saturated vacancies during the quench, do not normally contain stacking 
faults and are of dislocations with Burgers vectors a/2 [110]. The appear- 
ance of the copper sheets through which alpha-particles have passed (fig. 2) 
resembles that of a quenched foil. Closed loops can be identified but their 
shapes are less regular than those in quenched metals and in addition there 
is a very fine background of dots uniformly distributed throughout the 
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foil. The annealing behaviour is similar to that of a quenched metal 
(Silcox and Whelan 1960) the loops disappearing after several hours at 
350°c. The behaviour of the sheets containing the helium upon annealing 
is notably different and it is this which indicates the nature of the clustersy. 


3.1. Nature of Large Loops 

The larger dislocation loops (r~ 200A) which formed in the copper 
disappeared at about 350°c but at this temperature began to grow when 
helium was present. As the helium atoms are known to be very effective 
in capturing vacancies (Barnes et al. 1958, Barnes 1960) they will extract 
vacancies from any nearby dislocation by climb, if the temperature is 
sufficient for vacancy formation and migration. In copper this temperature 
is about 350°c. Dislocation loops formed by the clustermg of vacancies 
would contract in such a situation. Thus the growth of the large loops at 
this temperature suggests that they are formed by the clustering of inter- 
stitial atoms. The anneal also makes the loops more regular which would 
happen during climb. Their further growth to form a tangle of dislocation 
lines at 450°c is consistent with this. 

The total area of these large loops can be used to estimate the number of 
agglomerated interstitial atoms. The values given in the table, assuming 
that the loops lie on {111} planes, suggest that there are about 1-6 x 101* 
interstitial atoms cm? both in those thin films containing helium and those 
not. 

The annealing of the loops in regions containing no helium was such that 
the loops became smaller and disappeared after 2 hours at 350°c. Pre- 
sumably the vacancies which were present in equal numbers to the inter- 
stitial atoms, annihilated these loops. 

Near grain boundaries (fig. 26) there is a zone (~8x10-°cm wide) 
which contains none of the larger loops. Presumably the interstitial atoms 
have been captured by the boundary instead of forming loops, suggesting 


that they can normally migrate distances of this magnitude during the 
bombardment. 


3.2. Nature of the Dots 


The background of small dots, seen in fig. 2 (b), continues right up to the 
grain boundary, with no depleted zone, suggesting that they have formed 
from point defects which can diffuse only > 500 A during the bombardment. 
They must be very different from those producing the large loops and are 
probably formed from the vacancies which must be present in nearly equal 
numbers. 

It is very difficult to estimate the number of point defects clustered in 
these small dots which do not appear to be dislocation loops; their size 
is uncertain and may be determined by a strain field, also their number may 
pa en SE a 


} A difference in annealing of hardness between regions through which 
alpha-particles had passed and the region in which they were deposited had 
been noticed previously (Ghosh et al. 1960). 
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be underestimated because many are too small to be resolved. An estimate 
(table) suggests that if they are dislocation loops there are only about one 
fifth as many vacancies in them as there are interstitial atoms agglomerated 
in the larger loops. However, if the same values are used, but it is assumed 
that they are spherical voids, the number of vacancies agglomerated in 
them (table) more nearly equals the number of interstitial atoms, as should 
be the case. Certainly very small vacancy clusters would have a lower 
energy if arranged as a sphere rather than as a plate (Kuhlman-Wilsdorf 
and Wilsdorf 1960)*. 

The behaviour of the aluminium foils was similar to that of copper, 
but presumably in the former the vacancies and interstitial atoms had 
annihilated during the bombardment except where helium was present, 
where the appearance of the dislocation tangle and the formation of 
bubbles was almost the same as in the copper. 


3.3. Helium Bubble Nuclei 


After annealing at 450°C small gas bubbles can be seen as white dots 
near the grain boundaries on the photomicrographs (fig. 5b). Eventually 
after an anneal of | hour at 750°c these bubbles can be seen throughout the 
film if it is so oriented in the microscope that the dislocation lines do not 
give a contrast which obscures the bubbles (fig. 5d). These spots are 
40 A in radius and have the characteristics of gas bubbles ; they are spherical, 
i.e. they appear round with light centres becoming darker at their edges, 
with no diffraction contrast. The bubbles appear in rows and lie upon the 
dislocation lines (Barnes 1959-1960). This can be demonstrated by 
re-orientating the foil to give contrast effects at the dislocations. 

The number of helium bubbles observed once there are sufficient vacancies 
available is ~10''cm~* almost the same as the number of small dots 
(table). It seems possible that these dots are the nuclei upon which th 
helium bubbles form. 

The number of helium atoms (7) contained in the films can be calculated 
from the photomicrographs of the bubbles, using the formula (Barnes 1959), 


(2y/r) §rren = mkT 


where y is the surface energy (about 10% ergscm~*), r is the bubble radius 
(404), m is their number in unit area of the film (1-4 10"'cm~), & is 
Boltzmann’s constant, and 7' the absolute temperature (1023°K), therefore 
m= 1-4 x 104%cm-? (i.e. an atomic concentration of approximately 1-7 x 10-* 
in the copper). One concludes that there are approximately nine helium 
atoms for each interstitial atom in the loops and for each vacancy in the 


+ Recent observation of foils annealed for 1 hour at 400°c show that in some 
grains the remaining dots become triangular. This implies that the clusters are 
tetrahedra of stacking fault (Silecox and Hirsch 1959) both because of their shape 
and their stability above the self-diffusion temperature. 
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black dots. The small number of interstitial atoms observed in all 
these films (table) implies that most of the helium atoms were interstitial 
during the bombardment, because if each helium atom was substitutional 
1-4 x 10 vacancies em? must be provided and 1-4 x 10!° interstitial atoms 
cm-2 would be left ; only 1/10th this number was observed. 

The clusters would capture helium atoms during irradiation and gas 
bubbles would form upon them. However, presumably because of the 
shortage of vacancies these bubbles do not grow and thus cannot accommo- 
date further helium atoms. It can be estimated from the gas laws that 
only a fraction of the helium atoms reside in the vacancy clusters, during 
bombardment. Subsequent heating permits the acquisition of vacancies 
from the interstitial atom agglomerates, which grow. The number of 
vacancies these loops (radius ~ 200 A) produce during their climb to form 
the tangle of dislocations (equivalent to a radius of about 10004) is such 
that approximately two vacancies are provided per helium atom. This 
would be sufficient to enable each helium atom to become substitutional. 
The process ceases either because the helium bubbles eventually lock the 
dislocations and prevent their further climb, or because the dislocations 
become locked by mutual interaction. Whichever the reason, the structure 
remains stable until vacancies are provided from the grain boundaries 
(fig. 5e) and then the helium atoms produce observable bubbles as des- 
cribed previously (Barnes 1960). The dislocation tangle also enlarges 
with further heating and continues to enlarge as the bubbles grow larger 
and reduce in number. 


§ 4. CONCLUSIONS 

Copper bombarded with 1-4 x 10! alpha-particles em~? at about 200°c 
contains both dislocation loops r ~ 200 A formed by the clustering of inter- 
stitial atoms, and much smaller and more numerous dots (r~ 204) 
which are probably clusters formed from the less mobile vacancies. 
On annealing at about 350°C the clusters annihilate one another by the 
generation and diffusion of single vacancies so that the material reverts 
much to its original state. However, when the material contains large 
numbers of helium atoms these stabilize the vacancy clusters, and at about 
350°C they are able to grow by taking vacancies from the dislocation loops 
which, because they are formed from interstitial atoms grow, eventually 
becoming a tangle. The steady acquisition of vacancies by the small 
bubbles enables them to accommodate more helium atoms until eventually 
they can be clearly resolved in the electron microscope. 

Thus it appears that the vacancies produced by irradiation provide some 
space for the helium atoms to form bubbles. However, more vacancies are 
supplied from the dislocation loops which climb, but the total number is 
limited as Suggested previously (Barnes 1960) so that the size of the bubbles 
which form is too small to be detected using the normal optical microscope 


until large numbers of vacancies have been supplied from the grain 
boundaries. 
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ABSTRACT 


Recent observations of the growth of quenched-in dislocation loops 
during annealing are discussed. A mechanism for the coalescence of loops 
is proposed. The mechanism depends on short-circuit diffusion around the 
periphery of prismatic dislocation loops. This mechanism is contrasted 
with the theory of loop coalescence by vacancy diffusion through the matrix. 
It is shown that the experimental observations are in better accord with the 
short-circuit diffusion theory. 

Interpretation of the experimental observations on the basis of the 
short-circuit mechanism leads to a value for the activation energy for 
diffusion along a dislocation line. For aluminium this value is calculated 
to be 0-7 ev. 


§ 1. InTRODUCTION 


In recent electron-micrographic experiments on the annealing of quenched 
aluminium, Silcox and Whelan (1960) and Vandervoort and Washburn 
(1960) have observed that the density of dislocation loops present in the 
as-quenched material decreased markedly upon annealing. In the experi- 
ments of Silcox and Whelan, annealing was carried out in two distinct ways: 
annealing of the foil in situ within the microscope so that the details could 
be observed as annealing proceeded, and annealing of the bulk material 
before thinning, in which case the progress of annealing could not be 
followed continuously. 

It was found that the annealing of loops in foil became noticeable at a 
temperature of 170°c and proceeded by the shrinking of all loops, and that 
there was no tendency for any of the loops to grow at the expense of the 
others. Silcox and Whelan have accounted for these observations by 
proposing that vacancies are evaporated from the periphery of each loop 
and that these vacancies then diffuse to the surface of the foil. Since the 
free surfaces of the foil provide a nearby sink for vacancies, this process can 
occur rapidly at moderate temperatures. 

In specimens prepared from material annealed in bulk, however, they 
found that the decrease in loop density was accompanied by the growth of 
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some loops. This agrees with the observations of Vandervoort and 
Washburn, who found that considerable loop growth occurred in samples 
prepared from bulk material which had been annealed for ten minutes at 
150°c. Insamples which had been annealed for ten minutes at 125°c little 
loop growth was observed, but after ten minutes at 200°c annealing of the 
loops seemed to be essentially complete, with only a very few large loops 
remaining. 

That the growth of the loops must involve vacancy diffusion follows from 
the generally accepted view that the loops have originated by the collapse 
of vacancy precipitates. Accordingly, Silcox and Whelan have interpreted 
their observations on material annealed in bulk form in terms of the evapora- 
tion of vacancies from the periphery of a small loop, and the subsequent 
diffusion of these vacancies through the matrix to precipitate on a larger 
loop. The purpose of this paper is threefold: to propose a mechanism by 
which loop growth can occur without employing diffusion through the 
matrix between two loops, to contrast this process with the vacancy 
diffusion process proposed by Silcox and Whelan, and to interpret the 
experimental results in terms of both processes. 

The diffusion of vacancies through the matrix (bulk self-diffusion) is not 
the only means by which loops could coalesce. It is proposed that an effective 
mechanism for the coalescence of loops involves a diffusion short-circuit 
process. This is illustrated in fig. 1, in which atoms diffuse rapidly 
along the periphery of the loop. The rate at which diffusion can take place 
along this short-circuiting path is determined by the activation energy for 
diffusion along a dislocation line. The activation energy for diffusion 
along dislocation lines has been determined for various metals (Hoffman 
and Turnbull 1951, Hendrickson and Machlin 1954, Turnbull and Hoffman 
1954); the values obtained range between one-half and two-thirds of the 
activation energy for bulk diffusion. The activation energy for diffusion 
along dislocation lines in aluminium has not been measured, but it would be 
expected that it is appreciably lower than the activation energy for bulk 
self-diffusion. 

The elastic interaction between two nearby loops provides a driving 
force tending to create a net flow of diffusing atoms from the closest to 
the more distant edges of both loops}. As shown in fig. 2, this short-circuit 
diffusion around the periphery of the loops caused the two loops to move 
together. For convenience this process will be called ‘self-climb’. How- 
ever, self-climb alone can only account for the coalescence of coplanar 
loops. Loops can, in principle, also move by prismatic glide along a 
cylindrical surface inclined to the plane of the loop and containing the 
Burgers vector of the loop. The combination of self-climb and prismatic 
SoS test Leis dy ety: a hae eel eee i allen 


{The elastic interactions providing this driving force will be smaller in the 
foil due to the fact that a thin film of material can relieve these stresses by 
buckling. A comprehensive treatment of the limitations on dislocation 


interactions imposed by the thinness of foils used in electr i i 
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glide would therefore permit the motion of loops in any direction. The 
limitations on prismatic glide of loops have been discussed by Kuhlmann- 
Wilsdorf (1958), and will also be discussed below. 

Examination of the micrographs published by Vandervoort and 
Washburn (1960) shows that, in specimens prepared from material annealed 
at 150°c, the loop density is considerably smaller than the loop density in 
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Formation and migration of defects involved in self-climb. 


as-quenched material and, moreover, that the loops which remain after 
annealing are often considerably larger than those found, on the average, in 
the as-quenched material. In particular, a number of pairs of loops are 
observed. The distance between the members of such a pair is commonly 
smaller than the diameter of either loop. This effect is particularly evident 
in fig. 7 of Vandervoort and Washburn’s paper. From these observations 
the following conclusions have been drawn. 
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1. The coalescence of loops does not seem to be due to vacancy diffusion 
through the matrix. It is extremely difficult, for example, to understand 
how loops which are separated by distances such as those observed in the 
as-quenched material could have coalesced by diffusion through the 
matrix while some loops which are very near one another remained intact. 

2. The self-climb mechanism can satisfactorily explain these observa- 
tions. If this process controls the behaviour of loops during annealing, 
loops remain substantially intact until they finally merge, and this is 
consistent with the microscopic evidence. 


Fig. 2 


‘ , 

Successive stages in the coalescence of two coplanar loops by self-climb. The 
filled circles represent diffusing atoms and the open circles diffusing 
vacancies. 


However, another point must be considered, namely, whether the closely 
coupled pairs observed were just about to coalesce, or whether there are 
cogent reasons why they could not merge completely. 

It is possible for a vacancy plate on the (111) plane in a face-centred cubic 
metal to collapse in two different ways, depending on whether or not the 
atom planes above and below the plane of the plate are sheared along a 
direction lying in the plane of the plate as they come together. Ifno shear 
takes place, then Burgers vector of the resultant dislocation loop will be 
(a/3) [111] and the loop will contain a stacking fault since the stacking 
sequence normal to the plane of the loop is ABCACABC.... No stacking 
faults were observed in the specimens examined by Vandervoort and 
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Washburn or by Silcox and Whelan. If, on the other hand, the atom 
planes above and below the plane of the vacancy plate are displaced with 
respect to one another as they come together by one of the shears (a/6)[112], 
(a/6)[211], (a/6)[121], then the stacking sequence will be the normal 
ABCABCA ... so that the loop will not contain a fault, and the slip vector 
of the resulting loop will be (a/2)[110], (a/2)[011], or (a/2)[101]. 

These slip vectors are translation vectors of the face-centred cubic lattice 
so that a dislocation loop having any of these three slip vectors can glide 
on a cylindrical surface inclined to the plane of the loop and containing (as 
the generators of the cylinder) the slip vector of the loop. However, in 


Fig. 3 
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Geometries for calculating rates of loop coalescence. 
(a) Geometry for coalescence of coplanar loops by self-climb. The 
areas I and II are the interiors of the loops. 
(6) Geometry for coalescence of loops by bulk diffusion. 
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general not all parts of such a loop will be able to slip easily on this cylindrical 
surface since the surface will not be entirely composed of close-packed 
{111} planes. Such a situation is shown in fig. 4(a), in which the [110] 
segments of the loops can only slip on (001), while the segments oriented 
along [101] and [011] slip readily on the close-packed {111} planes. Since 
hexagonal loops corresponding to the configuration of fig. 4 (a) have often 
been found in foils, it is of interest to examine the possibility of slip on (001). 
Boas and Schmid (1931) have observed that macroscopic slip can occur 
on {100} planes in aluminium at temperatures above 400°c. It is possible 
that the attractive force between the neighbouring [110] segments might 
be high enough to cause them to combine by slip on (001) even at tempera- 
tures as low as 150°c. 
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For the geometry of fig. 4(a), where both loops have the same Burgers 
vector, the stress between the neighbouring [110] segments of the loops is 


peal AN ( 1 ) 
~ 2Qr(1—v) a/ 2d 


where pis the shear modulus, vis Poisson’s ratio, |b |is the magnitude of the 
slip vector, and d is the distance between the respective (111) planes on 
which the loopslie. To estimate the magnitude of this stress it is necessary 
to make some estimate of the value of d. Bearing in mind that the foils 
used by Vandervoort and Washburn and Silcox and Whelan were of the 
order of 10004 thick it seems reasonable to choose d~2004. Then 


Fig. 4 
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Coalescence of loops on parallel (111) planes by prismatic glide. 


(a) Combination of non-coplanar hexagonal loops. 
(b) Combination of loops capable of easy prismatic glide. 


o~5 x 10%dynes/cm?; this is roughly one hundred times as great as the 
macroscopic shear stress at which glide begins on {111} in aluminium. 
Since Boas and Schmid (1931) found {100} slip at 400°c under conditions 
for which the shear stress on planes of the {100} type was only slightly 
larger than the shear stress acting on planes of the {11 1} type, the possibility 
of inducing {100} slip at 150°c by a stress one hundred times higher than 
that required for {111} slip cannot be discounted. 

Prismatic glide of hexagonal dislocation loops has been observed by 
Hirsch et al. (1958). The observations were made on thin foils of quenched 
aluminium at temperatures lower than 100°c and provide strong experi- 
mental evidence for glide on {100} planes at ordinary temperatures. 
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As shown in fig. 4(b), it is possible to construct a loop all parts of which 
are oriented for easy slip on close-packed planes. Non-coplanar loops of 
this type can readily coalesce by self-climb followed by prismatic glide on 
{111} planes. 


Fig. 5 
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Successive stages in the coalescence of coplanar loops on (111) with different 
Burgers vectors. The dislocation segment formed at the intersection 
of the loops sweeps over the shaded area in (0). 


Coplanar loops having different slip vectors combine as shown in fig. 5. 
A segment of total dislocation is formed by the reaction 


(a/2)[110] + (a/2)[101]—> (a/2)[01T] 


where the loops join. The glide plane of this segment is the plane of the 
loops themselves. The segment sweeps over the interior of one or the 
other of the two loops and combines with the dislocation bounding this 
loop according to the reaction 


(a/2)[011] + (a/2)[110] > (a/2)[ 107], 


so that the slip vector of the resulting large loop is the slip vector of the 
‘unswept’ small loop. 
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§ 2. Kinetics or Loop CoMBINATION 


In the following section the kinetics of loop coalescence by bulk diffusion 
and by self-climb will be examined and compared. Rates of coalescence 
for loops having the size and separation typical of as-quenched aluminium 
will be estimated, and the observations of Vandervoort and Washburn 
(1960) will be used to derive an approximate value of the activation energy 
for diffusion along a dislocation line in aluminium. 

Figure 1 shows that the création of an atom capable of diffusing rapidly 
in the core of the dislocation loop is accompanied by the creation of a 
vacancy on the periphery of the loop. Coalescence of neighbouring loops 
by self-climb requires the diffusion of atoms along the loop cores to the 
farther edges of the loops and the diffusion of vacancies along the loops to the 
closer edges; thus, the loops approach one another by the transport of 
material around each loop from the closer to the farther edge. Although a 
diffusing atom is effectively constrained to remain in the core of a loop by its 
strong repulsive interaction with the atoms surrounding the core, a 
vacancy can diffuse away from the periphery of the loop into the matrix. 
It is this possibility of evaporation of vacancies from a loop which leads to 
the coalescence of loops by bulk self-diffusion. 

There are, however, two factors which argue against the evaporation of 
vacancies from the periphery of the loop. The first is that a vacancy 
created on the periphery of a dislocation loop must acquire an additional 
increment of energy to escape from the stress field of the loop. The second 
is that the rate of migration of a vacancy along the periphery of a dislocation 
loop should be much more rapid than the rate of migration away from the 
loop. In the absence of experimental values for the energy binding a 
vacancy to the compression side of a dislocation, and the activation energy 
for migration along a dislocation line, it is not possible to calculate the 
probability that a vacancy will escape from the periphery of a dislocation 
loop before it can diffuse along the loop and add to an incomplete row of 
vacancies. However, since the energy barrier a vacancy must surmount in 
diffusing away from the loop is an appreciable fraction of the difference 
between the activation energies for bulk diffusion and diffusion along a 
dislocation line, the probability of evaporation from the loop is smallat 
temperatures below about 150°c. 

The rate of transport of material from a source to a sink by means of 
diffusion is determined by the rate at which an individual carrier moves 
from source to sink and the number of carriers in the region in which diffusion 
is taking place. If it is assumed that the number of vacancies in the diffu- 
sion region is fixed (e.g. by requiring that thermodynamic equilibrium be 
maintained) then the number of vacancies diffusing from source to sink in 
time 7 is 

n= N(t/t) TE oe al 


where NV is the number of vacancies in the diffusion region and ¢ is the 
average time required for a vacancy to diffuse from source to sink. 
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In calculating the rate of coalescence of two neighbouring loops by 
either the bulk self-diffusion process or the self-climb process, it will be 
assumed that the initial concentration of vacancies is the concentration 
present in thermodynamic equilibrium at the annealing temperature, and 


that this concentration is maintained as the loops coalesce. 
If 


Us= the formation energy of a vacancy in perfect crystal, 


Um=the activation energy for migration of a vacancy in perfect 
crystal, 


U;' =the formation energy of a vacancy-atom defect pair capable of 
diffusing along a dislocation line, 


Uy’ =the activation energy for migration along a dislocation line, 


then U=U;+U ry is the activation energy for bulk self-diffusion and 
U’ =U?’ + Um’ is the activation energy for self-diffusion along a dislocation 
line. 

The time required for the coalescence of the two loops shown in fig. 3 will 
be calculated for the bulk diffusion process and for self-climb. For 
convenience one of the loops has been assumed to be much larger than the 
other, with the smaller loop having side a and the separation of the loops 
being J. 


§ 3. BuLK DIFFUSION 

On the basis of the bulk diffusion process it is to be expected that the 
smaller loop will shrink by evaporation of the (a/a))? vacancies from 
which it was formed, a, being the lattice constant of the crystal. These 
vacancies diffuse to the larger loop, causing to it grow. Referring to fig. 3, 
the distance these vacancies must diffuse is approximately 1. Assuming a 
random walk, the average time taken for a vacancy to diffuse from the 
smaller to the larger loop is therefore 


Ip Um 
oe Ree ae eua ty wide Pep Lai dt via (2 
L= Aine exp ‘at (2) 


where vy, is the vibrational frequency of the lattice. The volume in which 
diffusion is taking place is assumed to be approximately a*l, so that the 
number of vacancies in the diffusion region is 


al Us 
ne awe Poi irl a Ae (3 
IN he “exp ae 2 (3) 


Substituting ¢ from eqn. (2) and NV from eqn. (3) in eqn. (1) yields, for the 
time 7 required for the loops to coalesce, 


ra Nw 1 exp [Peel - exp tia} EGA 4) 


nt  Agvo Lovo kT’ 
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§ 4. SELF-CLIMB 


If the loops of fig. 3 coalesce by self-climb then both will move; however, 
due to the longer diffusion path around the larger of the two loops and the 
greater amount of material which must be transported around it, this loop 
will move much more slowly than the small one. Motion of the larger loop 
will therefore be neglected. 

The diffusion path around the periphery of the small loop is approxi- 
mately 2a in length. Assuming random walk and the same value vo for 
the atomic vibration frequency, the average time required for the transport 
of an atom from the near to the far side of the loop is 


4a? Bice 
isa . fyi ee SPIGA eS me ee 
t'~ Sex| = (5) 


In order for the loops to coalesce, all the material in region III of fig. 3 
must be transported around the loop so that if the quantities important 
in the self-climb process are denoted by affixing primes to the corre- 
sponding symbols in eqn. (1), then n’ ~al/a,?._ The number of atom-vacancy 
pairs present on the loop in thermodynamic equilibrium is 


4a U?' 
‘~ — -— 3 im opeiety eh actuate 
ar exp a (6) 
so that from (1), (5) and (6), the time required for the loops to coalesce is 
n't’ al Gey az { U' 
Ce ae OS ————— <p<i—}. ed 
Th N’ _ Uy? Vp exp { kT \ Ay? Vo exp a ( ) 


§ 5. COMPARISON wWitH EXPERIMENTAL RESULTS 


In the as-quenched material studied by Vandervoort and Washburn, 
the dislocation loops found were about 200—500 4 in size, and were separated 
by distances ranging from 1000 A to 50004. Bradshaw and Pearson (1957) 
have found U=1-20ev for aluminium; the value of U’ is not known. A 
reasonable estimate for the atomie vibration frequency is vy 10! /sec. 

Using the values a = 300 A, /= 3000 A, vp ~ 1018/sec, U = 1-20 ev, ay = 2°86 A, 
and various values of U’, the times r and 7’ have been calculated according 
to eqns. (4) and (7) for the temperature range of interest. These values are 
listed in the table. 

Since Vandervoort and Washburn found considerable loop growth 
following ten minute (6 x 10? sec) anneals at 150°o, inspection of the times 
listed in the table shows that : 


1. Bulk diffusion cannot account for the loop growth observed. 

2. Self-climb occurring with an activation energy of about 0-7ev does 
account for the observations. 

3. Since the time 7’ is much more sensitive to the value of U’ than it is to 
the geometrical factors in eqn. (7), the uncertainty in the value of U’ thus 
determined is relatively small. 


4. This value of U’ is consistent with the lowering of the activation 


energy for short-circuit diffusion processes that has been observed in other 
metals, 
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Approximate time required for combination of dislocation loops of 
diameter 300 A, separated by 3000 A 


Time required for combination of loops 


Temperature 
(°c) Bulk Self-climb 
diffusion = 
(U=1-20 ev) U’=0-6 ev 0-7 ev 0:8 ev 
Seconds 
100 levoe< 1K}! 125 2900 6-5 x 104 
125 1-6 x 10° 40 740 1-4 x 104 
150 2-0 x 104 14 220 3400 
175 3200 5:6 19 1000 
200 620 2°5 29 | 340 


§ 6. CONCLUSIONS 


1. The coalescence of quenched-in dislocation loops in aluminium takes 
place by short-circuit diffusion around the periphery of each loop. 

2. The activation energy for this short-circuit diffusion process is about 
0-7ev. (The measured value for the activation energy for diffusion 
through the aluminium lattice is 1-20 ev.) 
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ABSTRACT 


The magnetic susceptibilities of a number of solid solutions of chromium 
in vanadium have been measured at 20°, 77° and 293°K. The susceptibility 
at 20°K decreases with concentration from a value of 300 x 10-°e.m.u./g 
atom for pure vanadium to 159 x10-%e.m.u./g atom for pure chromium. 
Fine structure is observed in several regions of the susceptibility—composition 
curve. The temperature dependence of the susceptibility of all the alloys 
investigated was small. 

An interpretation of the susceptibility behaviour is given in terms of the 
collective electron model. 


§ 1. INTRODUCTION 


MAGNETIC susceptibility measurements can increase our understanding 
of the electronic band structure of metals and alloys. This applies 
especially to non-ferromagnetic transition metals—from which we have 
chosen vanadium—since a large part of their susceptibility is determined 
by the density of states of the conduction electrons at the Fermi level. 

The most important property of vanadium, however, which led to its 
selection for the present investigation is its ability to dissolve appreciable 
amounts of the majority of the transition elements in its own period. 
The concentration ranges (in atomic per cent) over which these elements 
form stable primary solutions with vanadium are (Rostoker 1958): 
titanium 0-70, chromium 0-100, manganese 0-50, iron 0-30, cobalt 
0-8, nickel 0-8 and copper 0-7. A comparison of the results for such 
vanadium alloys should, in principle, distinguish the effects of atomic 
concentration and electron transfer on the Fermi level and the shape 
of the conduction band. 

In this paper we report the results of magnetic susceptibility measure- 
ments for solid solutions of chromium in vanadium. Subsequently we 
will compare these results with those for the other alloy systems. 


§ 2. PREPARATION OF THE SPECIMENS 


The vanadium and chromium used for the alloys were obtained from 
a number of sources, purer materials being used during the course of the 
work as they became available. No systematic dependence of suscep- 
tibility upon the purities of these source materials was found for vanadium 


+ Communicated by the Authors. 
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or the vanadium-rich alloys containing up to about 30at. °%% chromium. 
This was not the case, however, for chromium and so the chromium-rich 
alloys, which were all prepared towards the end of the investigation, were 
made using only the purest materials. 

The suppliers of the vanadium were Johnson, Matthey and Company 
(J.M.), the Battelle Memorial Institute (B.M.L.) and the U.S. Bureau of 
Mines (U.S.B.M.). The methods of preparation and the purities (non- 
metallic elements being included in the analyses), were: 


J.M.: calcium reduction of vanadium pentoxide (99-63 weight %), 
U.S.B.M. : electrolytic refinement of vanadium metal (99-85 weight %), 
B.M.I.: Van Arkel process from vanadium iodide (99-92 weight %). 


Detailed analyses have been given in an earlier paper (Childs ef al. 1959 
referred to henceforth as [). 

The chromium was provided by Johnson, Matthey and Company 
(99-65 weight °%) and the Fulmer Research Institute (99-99 weight %). 
The latter material (F.R.I.) had been hydrogen-reduced to lower its 
oxygen content to 0-0032 weight %. 

Alloy buttons of about 3g weight were prepared by melting together 
appropriate amounts of vanadium and chromium in an argon-are furnace. 
The buttons, which were supported on a water-cooled copper hearth, 
were inverted and remelted a number of times so as to ensure adequate 
mixing of the two components. Contamination by oxygen or nitrogen 
during these melting operations was minimized by frequent gettering of 
the argon atmosphere with molten zirconium. A leak test on the 
apparatus was made before melting was begun. 

After melting, the buttons were annealed, both to reduce small-scale 
concentration fluctuations caused by segregation on solidification of 
chromium-rich material, and to promote the formation of a random solid 
solution. Each button, supported on a plate of purified thoria, was 
sealed in an argon-filled silica tube and heated by high-frequency 
induction for one hour at 1600°c. A leak test and argon gettering 
operation was performed before the annealing. 

The density and hardness (§3) of the button were then measured, 
after which susceptibility specimens approximately 8x 1-5x1-5mm in 
size were cut from its centre using a clean Carborundum slitting wheel. 
A stream of coolant liquid (a water-oil emulsion) prevented excessive 
heating and oxidation of the specimens during this process. The 


specimens were finally degreased in alcohol, cleaned in dilute hydrochloric 
acid and washed in distilled water, 


§ 3. EXAMINATION OF THE ALLOYS 


Before a specimen was accepted for susceptibility measurement it was 
carefully examined in order to determine whether its composition, 


homogeneity and purity were satisfactory. This involved the following 
procedure. 
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3.1. Composition 


The composition of each alloy button was determined from the initial 
weights of its constituents after correcting for the loss in weight which 
occurred on melting. The loss in weight, which was usually small, was 
assumed to be due almost entirely to chromium evaporation since the 
vapour pressure of chromium below 2000°c is more than one thousand 
times that of vanadium. 

These composition values were confirmed for certain selected specimens 
by direct chemical analysis. 


Fig. 1 
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Variation in density with composition for the vanadium-chromium system. 


alloys of J.M. vanadium and J.M. chromium; 
alloys of B.M.I. vanadium and J.M. chromium; 
alloys of B.M.I. vanadium and F.R.I. chromium, 
alloys of U.S.B.M. vanadium and F.R.I. chromium. 


Meope 


The final composition values are estimated to be correct to within 
either +0-5at.% chromium or +2% of the quoted chromium content, 
whichever is the lesser. 
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3.2. Density 


The density of each annealed button was measured at room temperature 
using the standard water displacement method. The variation in density 
with composition for the vanadium-chromium solutions is shown in 
fig. 1. Within the accuracy of measurement the density values for all 
the specimens fall on a common curve. 

The density values served as a further check on composition. They 
were also used (§ 6.3) in estimating the effects on the susceptibility values 
of changes in the interatomic spacing. 


Fig. 2 
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Hardness of the vanadium-chromium alloys as a function of composition 


alloys of J.M. vanadium and J.M. chromium: 

alloys of B.M.I. vanadium and J.M. chromium: 

aie Me nee vanadium and F.R.I. chromium; 

alloys of U.S.B.M. vanadium and F.R.I. i 

S.E.R.L., Baldock, vanadium. spit a 
Aeronautical Research Laboratories, Melbourne, chromium 
x---x_ Carlson and Eustice (1959). : 
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3.3. Hardness 

After measuring its density the button was cut in half and the cut 
surface of one piece polished. The Vickers hardness of this piece was 
then measured at a number of points across the surface. 

Any large-scale variation in concentration across the specimen (which 
was never found) would have produced a corresponding variation in 
hardness, whereas excessive contamination throughout the specimen 
(for example by oxygen) could be detected by the resulting high average 
hardness. 

Figure 2 shows the hardness of a number of specimens as a function of 
their chromium concentration. Also plotted in this figure are hardnesses — 
for vanadium, reported in I, and some additional values obtained from 
the results of Carlson and Eustice (1959). Their published results have 
been converted from Rockwell A to Vickers hardness values by Dr. O. N. 
Carlson (private communication). 

The hardness of alloys prepared from the purest material—B.M.I. or 
U.S.B.M. vanadium and F.R.I. chromium—increased from about 90 VPN 
for pure vanadium to 400VPN for 50% chromium. It then fell to 
150 VPN for pure chromium. The hardness values of the alloys prepared 
from the less pure materials occupied a band approximately 150 VPN 
wide above the base curve. 


3.4. Microstructure 


The polished and etched surfaces of the buttons were examined at 
high magnification (x 500) for traces of foreign material dispersed as a 
precipitate or in other form. Most of the specimens, in fact, contained 
such a dispersion usually as a fine filamentary network within the grains 
with, in some cases, isolated particles distributed at intervals along the 
grain boundaries. Typical micrographs (x 250) of the alloys are repro- 
duced in fig. 3 (Pl. 165). 

The precipitate is presumably associated with the impurities initially 
present in the alloy constituents since its volume was generally much less 
for alloys made from the purest materials. Experiments with aged 
specimens (see § 5) showed, at least for the vanadium-rich alloys, that the 
susceptibility was not significantly altered unless the precipitate occupied 
more than one volume percent of the specimen. Such specimens were 
therefore rejected. 

§ 4. SUSCEPTIBILITY MEASUREMENT 


The method used to measure the susceptibility of the specimens has 
been described previously (Childs and Penfold 1957). The measurements 
were made using a Sucksmith ring balance to an estimated accuracy of 
+1% relative to the susceptibility of tantalum (0-849 x 10~°e.mu./g) 
as standard. The apparatus was frequently re-calibrated during the 
course of the work. The measurements were made at three temperatures, 


20°, 77° and 293°K. 
482 
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§ 5. RESULTS 


The mass susceptibility values are given in full in the table. This 
table includes some results for a pure vanadium sample (S.H.R.L.) 
which were not reported in I. This specimen was produced by the Van 
Arkel process and further refined by zone-melting at S.E.R.L. Baldock. 
The exceptionally low hardness of the specimen (57 VPN) indicated 
that ite purity was greater than that of any other vanadium sample 
investigated. 


The mass susceptibilities of vanadium—chromium solid solutions 


Chromium Xm X 10% e.m.u./g 


concentration pene at 
at. % Vanadium Chromium 20°K 77°K 293°K 
0 J.M. — 5-94 5-96 5-83 
0 U.S.B.M. = 5°83 5-86 5:77 
0 B.M.I. — 5-89 5-88 5-81 
0 S.E.R.L. — 5-92 5:95 5-86 
1-00 J.M. J.M. 5:87 5-83 5-63 
1-57 B.M.I. J.M. 5-64 5-67 5-55 
1-91 B.M.I. F.R.1. 5-59 5-62 5-51 
1-96 J.M. J.M. 5:64 5-63 5-57 
1-96 U.S.B.M. F.R.1. 5-46 5-50 5:37 
2-84 J.M. J.M. 5-64 5:64 5-57 
3-90 U.S.B.M. F.R.I. 5-82 5-72 5-64 
4-00 J.M. J.M. 5-74 5-68 5-75 
6-9 J.M. J.M. 5-62 5-63 5-54 
10-5 J.M. J.M. 5-65 5-63 5:57 
12-2 J.M. J.M. 5:30 5:30 5-27 
15:3 US.B.M. F.R.I. 5-27 5:28 5-26 
19:2 J.M J.M. 5-49 5-45 5:43 
29-1 B.M.I. J.M. 5:05 5-04 5-02 
38-5 U.S.B.M. F.R.I. 4-84 4:79 4-79 
48:8 U.S.B.M. F.R.I 4-88 4-90 4-76 
50-0 U.S.B.M. F.R.I 4:79 4-81 4-73 
50-5 U.S.B.M. F.R.I 4-91 4-89 4-76 
59-5 U.S.B.M. F.R.I 4-04 4-04 4-04 
69-2 U.S.B.M. F.R.I 3°51 3-51 3-58 
77:8 U.S.B.M. F.R.I 3:27 317 3-23 
100-0 — J.M. 3:38 3-47 3:40 
100-0 — F.R.I 3-01 2-99 3°19 
100-0 — A.R.L 2:99 2-98 3°17 


The estimated absolute accuracy of Xm is +0-10 x 10-*§ e.m.u./g. 


The atomic susceptibilities at 20°K are shown as a function of the 
chromium concentration in fig. 4. Also included in this figure are some 


results obtained by Lingelbach (1958) at 100°x for alloys of high chromium 
concentration. 
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The suceptibility at 20°K decreases with concentration in an oscillatory 
manner up to 50at.% of chromium and then falls sharply. Thus there 
are small minima at 2 and l5at.%, with maxima at 6 and 20at. % 
chromium. In addition there is a broad minimum at 32 at. % and a 
relatively broad maximum at 50 at. % chromium. 


Fig. 4 
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Atomic susceptibility at 20°K for the vanadium-chromium system as a function 
of composition. 


alloys of J.M. vanadium and J.M. chromium; 
alloys of B.M.I. vanadium and J.M. chromium; 
alloys of B.M.I. vanadium and F.R.I. chromium; 
alloys of U.S.B.M. vanadium and F.R.I. chromium. 
Lingelbach (1958) 100°K. 


XMeop 


In order to confirm this fine structure several alloys were investigated 
having compositions of about 2at.% chromium and 50 at. % chromium. 
Four separate specimens having a nominal 2at.°%, were prepared from 
various combinations of the vanadium and chromium source materials 
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to determine whether this minimum was a genuine alloy effect unconnected 
with the presence of small amounts of impurites. The hardnesses of 
these specimens ranged from 110-314VPN yet, as shown in the figure, 
their susceptibility values are grouped closely together and show no 
systematic dependence on hardness or on the purity of the source 
material. This satisfactory agreement suggested that it was not 
unreasonable to continue the curve through all our observed points 
and so obtain the rest of the fine structure. 

The relatively broad maximum at 50 at. % chromium may be associated 
with the development of long-range order in this alloy. To test this 
possibility two specimens were prepared from buttons which were 
quenched rapidly from 1900°c to room temperature—the normal 1600°c 
annealing stage, followed by a relatively slow cooling of the alloy, being 
omitted. The susceptibility of these quenched specimens was the same, 
within experimental error, as that of the annealed and slowly cooled 
specimen. Thus ordering, if it occurs, must reach equilibrium in a 
period of less than a few seconds. 

Some further evidence on this point was provided by a neutron 
diffraction study of the annealed specimen. No significant superlattice 
line (100) associated with the ordered state was detected. From an 
estimate of the errors of measurement it was concluded that the amount 
of long-range order present must be less than 10°. This result is expected 
since the proximity of vanadium and chromium in the electrochemical 
series and their similar atomic sizes do not energetically favour ordering. 

Since the major emphasis in this work was on the behaviour of vanadium 
as a solvent metal no attempt was made to investigate the suscepti- 
bility variation for concentrations above 78 at. 94 chromium. However 
Lingelbach’s (1958) results, some recent resistivity data (§ 6) and electronic 
specific heat measurements of Wei, Cheng and Beck (private com- 
munication) show that further fine structure in this region of the 
susceptibility curve may be expected. 

The variation in susceptibility (table) between 20° and 293°xK for all 
the alloys was very small and not appreciably dependent on concen- 
tration. On the average the susceptibility decreased by 1% over this 
temperature range. The absence of any strong temperature dependence 
for the 78at.% chromium alloy is noteworthy since, for an S0at. ve 
chromium specimen, Lingelbach found the susceptibility to drop from 
227 x 10% e.m.u./g atom at 66°K to 180 x 10-%e.m.u./g atom at 293°K—a 
20% variation. In view of the lack, in our results, of any systematic 
trend in the temperature coefficient of susceptibility between 0 and 185 
chromium, it seems likely that this particular specimen of Lingelbach’s 
was impure. 

As mentioned earlier, the susceptibility of the J.M. chromium differed 
appreciably from that of the F.R.1. material, the values at 293°K being 
177. and 166 x 10e.m.u./g atom respectively. A third chromium 
sample, obtained from the Aeronautical Research Laboratories, Melbourne, 
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which also contained only 0-003 weight °%% oxygen, had a susceptibility 
(165 x 10-%e.m.u./g atom at 293°K) close to that of the F.R.I. sample. 
The higher value for the J.M. specimen is presumably due to the 0-37 
weight % of impurity—mostly oxygen—contained in it. 

In the case of vanadium and of the alloys with up to 30 at. % chromium, 
excessive contamination by oxygen and nitrogen (introduced either 
deliberately or by accidental mishandling during preparation) caused 
the susceptibility to decrease. This behaviour, which is opposite to that 
shown by chromium, has been mentioned in I for a sample of pure 
vanadium. ‘This specimen decreased in susceptibility from 296 to 
244x 10% e.m.u./g atom on solution of 7 weight % of oxygen. 

No significant difference was found between the susceptibility values 
for specimens prepared from the different source materials and having 
hardness values between the two curves shown in fig. 2. Samples with 
hardnesses much above the upper curve gave anomalously low suscep- 
tibility values and were rejected. 

Some specimens which had, after the initial measurements, been 
allowed to stand at room temperature for more than one year showed an 
ageing effect whereby the susceptibility decreased with time. One such 
specimen—the 10-5at.°% chromium sample—decreased in susceptibility 
from 284 to 273 x 10-*e.m.u./g atom at 293°K after 28 months. Removal, 
by etching, of a surface layer 1:-4x10-%cm thick did not affect the 
susceptibility of the specimen, thus indicating that the ageing process 
is a volume rather than a surface effect. This was confirmed by cutting 
a new specimen from the middle of one of the halves of the original alloy 
button. This new specimen had, within experimental error, the same 
diminished susceptibility. 

This effect is presumably connected with impurity precipitation since, 
while the 10-5at.°% alloy originally appeared on metallographic 
examination perfectly free from inclusions—the grain boundaries being 
sharp and fine—re-examination after ageing revealed the presence of a 
substantial quantity of precipitate distributed preferentially along the 
grain boundaries but also throughout the grains. 


§ 6. Discussion 


Before attempting to explain the variation with composition of the 
alloy susceptibilities, it is necessary to consider which of many conceivable 
magnetic processes contribute appreciably to the susceptibilities of the 
pure component metals. 

6.1. Vanadium Metal 


For vanadium the major contribution to the susceptibility most 
probably. comes from the spin paramagnetism of the outer electrons (five 
per atom). If any of these electrons were in localized atomic-like orbitals 
then, from Hund’s rule, their spins would necessarily be parallel. In 
‘consequence each atom would show a permanent magnetic moment, 
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However, the results of neutron diffraction measurements made at 20°K 
by Shull and Wilkinson (1953) show that, within the limits of experimental 
error (0:1 Bohr magneton), there are no such moments, aligned or 
unaligned, on the vanadium atoms. All the outer electrons must there- 
fore be in non-localized states in the conduction band of the metal where 
they will give rise to a Pauli-type paramagnetism. 

Conduction electrons, if mutally non-interacting, contribute a term 
to the susceptibility per gram atom (Pauli 1927) given by 

Xp = 2B’ Nn(Ey) 
where N is Avogadro’s number and £ the Bohr magneton. n(H)), the 
density of states for one direction of spin per unit energy range per atom 
at the Fermi surface, may also be expressed by the relation (Mott and 
Jones 1936) 

3y 

to) = 5 81EN 
where k is Boltzmann’s constant and y is the temperature coefficient of 
the electronic specific heat per gram atom. Substituting for y the value 
8-8x103J/g atom deg? obtained by Corak et al. (1954), we find that 
for vanadium n(H,)=1-86 states/ev atom. The Pauli susceptibility for 
vanadium is thus 122x10-*e.m.u./g atom. The observed value (x9) 
at 20°K, corrected (see later) for diamagnetic terms is 308 x 10-*e.m.u./g 
atom. 

Such a large difference between y, and y, is not peculiar to vanadium 
but is observed in all transition metals. It has been attributed (e.g. 
Mott 1935) to a spin-interaction of a ferromagnetic nature between the 
conduction electrons. This ‘exchange interaction’ tends to increase 
the unbalance produced by the field in the number of parallel and anti- 
parallel electron spins and hence to increase the induced moment. 

If we write (Stoner 1936) 

1 ade 

Xo Xp 
then K, which is proportional to the exchange interaction, has a value for 
vanadium of 5000¢atom/e.m.u. 

There are two small diamagnetic susceptibility contributions to be 
considered. The first of these, which is associated with the orbital 
motion of the tightly bound electrons on the vanadium ions, has a value 
(Van Vleck 1932) of 7-7 x 10-%e.m.u./g atom, i.e. 2-5°% of the susceptibility 
at room temperature. The second term arises from the curvature imposed 
by the magnetic field on the paths of the conduction electrons. This 
term cannot be calculated exactly but is probably negligibly small as the 
following rough argument shows. 

The value for y appropriate to a free electron gas containing five elec- 
trons per atom is 0-94x 10-J/g atom deg? (Mott and Jones 1936). 
Vanadium metal, as mentioned however, has a y of 8-8 x 103 J/g atom 
deg” so that the effective thermal mass, mt, of the electrons at the Fermi 
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surface is about nine times the normal electronic mass, m. We shall 
take m+ equal to ma, the effective mass of the electrons accelerated by the 
magnetic field (Cohen 1958). Then, since the conduction electron 
diamagnetism (Landau 1930) is associated with the electronic orbits in 
momentum space it is inversely proportional to m2 and not m2. Thus, 
instead of being one-third of the paramagnetic susceptibility, as is the 
case for free electrons, this diamagnetic component for vanadium is only 
0:4% of the paramagnetism. 


6.2. Chromium Metal 


Chromium is much more complex than vanadium in its magnetic 
behaviour. Shull and Wilkinson (1953) and later Corliss et al. (1959) 
found the metal to be antiferromagnetic with a mean moment per atom 
of 0-4+0-058. The magnetic susceptibility of an oxygen and nitgrogen- 
free polycrystalline specimen was measured over the temperature range 
- 66° to 805°K by Lingelbach (1958). His graph shows a marked kink at 
a temperature (310°K) significantly close to the Néel temperature (7'y) 
of 308°K observed by Corliss et al. for a single crystal. 

However, these results disagree with the value of 7’, observed by 
Shull and Wilkinson, namely 475°K, on a polycrystalline specimen. This 
suggests that 7’, may well depend on either the preparation or the 
impurity content of the sample. However it is difficult to reconcile this 
result with the fact that several of the properties of chromium show an 
anomaly only at about 310°K (Sully 1954). In particular, the heat 
capacity of a polycrystalline sample shows a lambda-type anomaly at 
311-5°K (Beaumont et al. 1960). We shall therefore assume in our sub- 
sequent discussion of chromium that 7'y is 310°K. 

In contrast with the usual behaviour of antiferromagnetic materials 
the susceptibility above 7’, does not decrease according to a Curie-Weiss 
law. Instead a progressive increase occurs, the value changing from 
166 x 10-*e.m.u./g atom at 310°K to 177 x 10-%e.m.u./g atom at 800°K. 

The electronic specific heat measurements made by Estermann et al. 
(1952) show that y for pure chromium is 1-54 x 10-°J/g atom deg”. 
The density of states at the Fermi level calculated from this value is 
0-33 states/ev atom—considerably less than that for vanadium. Wei 
et al. (1959) have found that y increases rapidly as manganese and iron 
are dissolved in chromium. It seems likely therefore that the Fermi level 
for this metal falls just above a minimum in the density of states curve. 
This suggestion is well-supported by the measurements of soft x-ray 
absorption in thin films of chromium by Agarwal and Givens (1957). 

The Pauli susceptibility of chromium at 0°K calculated from the above 
density of states value is 21x10-%e.m.u./gatom. If we assume that 
the mean moment of each chromium atom above the Née! point is 
0-48 and the Curie-Weiss 9, like T'y, is 310°k, then the corresponding 
Curie-Weiss susceptibility at the Néel point is 194 x 10~%e.m.u./g atom. 
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The sum of the Curie-Weiss and Pauli terms at 310°xK (215 x 10~°e.m.u./g 
atom) does not agree with the measured susceptibility value (173 x 1076 
e.m.u./g atom), corrected for the ion core diamagnetism (6-6 x 10-*e.m.u./g 
atom—Van Vleck 1932). This is not particularly significant since a slight 
increase in the value of @ could produce agreement. However, no 
change in @ or the mean moment per atom can account for an increase 
in the susceptibility of chromium above the Néel temperature. Thus 
according to the above figures the susceptibility should decrease by 
87x 10-%e.m.u./gatom between 310° and 800°K instead of increasing 
by 11 x 10-%e.m.u./g atom as observed. 

It thus seems possible that, above the Néel temperature, there are no 
spins in chromium obeying a simple Curie Weiss law and that a collec- 
tive electron model, such as the one developed by Lidiard (1954), is more 
appropriate. 

Assuming this to be the case, and extrapolating to 0°K the suscep- 
tibility values observed for temperatures above 7'y, then the corrected 
value of the susceptibility (y,) in the absence of antiferromagnetism is ~ 
170 x 10-e.m.u./gatom. Taking x, as 21 x 10-%e.m.u./gatom the value 
obtained for K is 42000gatom/em.u. The exchange interaction 
parameter is thus about eight times as large for chromium as for vanadium. 

The calculated K values for 12 other transition metals for which the 
appropriate experimental data are available are less than 20000g 
atom/e.m.u., except for tungsten (24000) and molybdenum (37000). 
The fact that the interaction parameter is large for these other two group 
VI metals, lends support to the idea that it should also be large for 
chromium. 

As we have implied the two diamagnetic contributions to the 
susceptibility of chromium are, as for vanadium, both small. Thus the 
effective thermal mass of the electrons at the Fermi surface in chromium, 
calculated in a similar way to that. described for vanadium, is about 
1-6 times the electronic mass. The Landau diamagnetic component is 
therefore about 13% of the Pauli paramagnetism and less than 2% of 
the observed paramagnetism at room temperature. The core dia- 
magnetism is about 4% of this observed figure. 


6.3. Vanadium=Chromium Solid Solutions 


The fall in susceptibility with composition for the alloys up to 80 at. % 
chromium is most simply explained in terms of a collective electron model. 
Thus the fact that the susceptibilities of all the alloys, like that of 
vanadium, are only weakly temperature dependent is evidence that no 
free permanent moments appear on either the vanadium or the chromium 
atoms. ‘This is supported by some resistivity measurements (De Vries 
1959) which show that the Néel temperature of chromium is very 
dependent on small alloying additions. De Vries found, for example, 
that the addition of lat. of vanadium decreases Ty by 85°K, which 
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suggests that there should be no moments present on the chromium 
atoms in an 80at.°% chromium alloy. 

If the structure of the vanadium conduction band is assumed to remain 
essentially unchanged on adding chromium, then the major effect of 
alloying will be to produce an increase in the Fermi energy as electrons 
from the chromium atoms enter the common band. The Fermi level 
for vanadium comes at a point where the density of states is high but 
decreasing (Kriessman and Callen 1954) while, as mentioned, that of 
chromium probably occurs just above a minimum in the band. Thus, 
on alloying, the Pauli paramagnetism should initially decrease and then, 
finally, increase slightly. 

On the basis of this model the observed susceptibility curve would 
represent the density of states curve as modified by the exchange 
interaction effect. The exchange parameter also will vary with 
composition. A possible relationship involving K is suggested by its 
behaviour, mentioned previously, for a number of transition metals. 
K for these metals is roughly inversely proportional to y,. If this were 
also true for these alloys then the observed susceptibility—concentration 
curve would be a direct representation of the density of states curve. 
All the observed oscillations would then be reflections of fine structure 
in this curve. 

Such oscillations would not produce a temperature variation in 
susceptibility greater than that observed between 20 and 300°K. In 
this temperature range the edge of the Fermi distribution is sufficiently 
smeared out so that only an average of the sharp oscillation at 2at. %, 
for instance, will be observed. However, by assuming that K is 
temperature independent a simple calculation shows that below 20°K 
the susceptibility of the 2at. % alloy may further decrease by about 3%. 

A small correction to the above model is necessary to allow for the 
increase in width of the band which must occur as a result of the decrease 
in interatomic spacing on alloying (§3.2). For example, if the lattice 
parameter of vanadium were reduced from 3-03A to 2-884 (the value 
for chromium) then the Pauli susceptibility, assuming a parabolic density 
of states curve, would decrease by 10%. ‘The effect accounts for about 
one tenth of the observed susceptibility decrease between vanadium and 
chromium. 

A confirmation of some of these ideas must await measurements of the 
electronic specific heat of all these alloys. It is gratifying that the results 
obtained at a few concentrations for the vanadium—chromium system 
by Wei, Cheng and Beck (private communication) qualitatively support 
our conclusions. 
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ABSTRACT 


The hyperfine structure (H.F.S.) of the J=3/2 ground state of 2°°Bi has 
been found by the atomic beam magnetic resonance method. Two of the 
three H.F.S. intervals were measured directly, giving 


(W50— W 60) /h=2884-7 +0-2 Me/s 


and (W4,—W5; 5)/a=2171-5+0-1 Me/s. These results are consistent with a 
magnetic dipole interaction constant A of —446-97 + 0:04 Me/s and an electric 
quadruple interaction constant B of —303:3+0°3 Me/s. The latter yields an 
uncorrected nuclear electric quadrupole moment of —0-34 barns, and the field 
dependence of A’ =0 transitions gives — 1-6433 + 0-0002 for the ground state 
gz; value. The results are discussed in terms of the intermediate coupling 
theory of Breit and Wills (1933). 


§ 1. INTRODUCTION 


Goudsmith and Back (1927) studied the H.F-.S. of the ?°°Bi ground state by 
optical means and showed that the ground state has a total width of about 
0:08cm-, or 2400 Mc/s, but they were unable to resolve the structure. 
Later Mrozowski (1942) resolved the largest of the three H.F.S. intervals 
and estimated a total width of about 6750 Mc/s; he also showed that the 
H.E:S. is inverted, i.e. the energy level with the smallest value of / has the 
largest energy and A is negative (see fig. 1). We have measured two of 
these intervals directly and the results are consistent with a total width of 
6641-2+0-7Mc/s. Rochester (1960) has confirmed that the sign of A is 
negative by the method described in the following paper on ™1Sb. 

The principles of the method, which was first applied to atoms by 
Kellog e¢ al. (1939), are given in detail in the book by Ramsey (1956). 
The atomic beam passes in turn through an inhomogeneous deflecting 
field H,, a homogeneous field H, and a second inhomogeneous deflecting 
field H, before falling on the detector (fig. 2). A weak oscillating magnetic 
field H, superimposed on the steady field H, induces transitions between 
certain states W, and W, with greatest probability when hy=|W,-— W,]. 
The deflecting fields are arranged so that only atoms which undergo 
suitable transitions in the C field, the ‘flopped’ atoms, are refocused on 
the detector, the others being deflected away from the axis. An increase 
in detected beam then indicates the resonance condition and thus allows 
measurement of the energy interval. 


+ Communicated by the Authors. 
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§ 2. ENercy LEVELS 


The energy levels expected for 2°Bi (J=3/2, [=9/2) in a steady 
magnetic field H, are shown schematically in fig. 1, which applies when the 
H.F.S. is inverted. The general way in which the energy of any level 
W(F, m,) varies with H, and depends on the constants J, J, f, mz, g;, A 
and B can be calculated by methods described by Ramsey (1956). Per- 
turbation theory, for example, yields the following expressions for fre- 
quencies associated with the AF =0 transitions marked by arrows in fig. 1: 
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The hyperfine structure of 2° Bi in weak and strong fields (schematic). 
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where #=9,u)H,/h and y=g,u9H,/h. The interaction constants A and B 
are related to the zero field intervals by 
vy =6A + 2B/3, 
vo = 5A —5.B/23,- Me/s. oD ee eat (3) 
vg=4A —2B/3, 
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Schematic arrangement of the atomic beam magnetic resonance apparatus. 
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§ 3. APPARATUS 


The apparatus used by Bellamy and Smith (1953) for radioactive alkalis 
was modified so that an atomic bismuth beam could be produced and 
detected. The ‘flop-in’ magnet system was arranged to refocus on the 
detector only those atoms in which the oscillating field caused transitions 
such that the sign of m, was changed and the selection rules for magnetic 
dipole transitions AF=0, Am,= +1 or AF=+1, Amp=0 or +1 were 
satisfied. The transitions marked in fig. 1 are suitable and these were 
studied in some detail. 

The atomic bismuth beam was produced in the oven shown in fig. 3, 
which is similar to that used by Leu and Fraser (1928). Electron bom- 
bardment of the molybdenum snout by a current of 100ma at 1400V 
prvduced a snout temperature of 200°K, sufficient to cause almost 
complete dissociation of the molecules coming from the oven cavity, 
which remained at a temperature nearer 1000°K. A 0-004in. wide oven 
slit was machined in the end of the snout. 

The beam was detected with an electron bombardment detector (Wessel 
and Lew 1953). Positive ions produced by an electron current of 4mA 
were analysed and counted by the mass-spectrometer and electron- 

‘multiplier used earlier by Bellamy and Smith (1953). The count-rate 
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increase on resonance was only about 3% because of a large ionizer back- 
ground so it was advantageous to amplitude modulate the beam or the 
radio-frequency field and use a phase-sensitive detector. Total beam 
intensities normally used were estimated to be about 10’ atoms per second 
at the detector. 

The homogeneous field H, was calibrated as required by observing 
resonances with a potassium beam from a separate resistance heated oven 
mounted beside the bismuth oven, and this beam was detected by surface 
ionization. Data published by Kusch and Taub (1949) were used to 
calculate the magnitude of H, for a given *°K resonance frequency. 


Fig. 3 


“ 


4 
1/4 


The bismuth oven. 


The oscillating field was produced by a simple slotted hairpin in the C 
magnet gap. It could be positioned to stimulate o transitions (Am,=0, 
H, parallel to H,) or z transitions (Am,= +1, H, perpendicular to H,) as 
required. Typical resonances were about 0-15 Mc/s wide and they were 
not broadened significantly by C field inhomogeneity. The o resonances 
were double peaked, as one expects with this type of loop (Woodgate and 
Hellwarth 1956). Normal heterodyne methods were used to set the 
oscillator frequencies to about 2 parts in 10°, 


§ 4. RESULTS 


The AF'=0 transitions v, and v, indicated in fig. 1 were observed at each 
of the magnetic fields shown in the table. Insertion of these frequencies 
and fields into eqns. (1) and (2) then gave the values of y, and v, shown in 
the table, where it is seen that their accuracy increased rapidly with the 
magnitude of H,. The nuclear magnetic resonance value of + 4:0797 NM 
obtained by Protor and Yu (1950) was used for ply: 

Having determined the H.F.S. intervals to about + 20 Me/s it was then 
possible to look for AF'= +1 transitions and so determine the intervals 
with greater accuracy. Microwave triode oscillators were used to provide 
pl power in the frequency range required. Fifteen z transitions and nine 
o transitions satisfy both the selection rules and the refocusing conditions 
in the apparatus, and all of these were observed in a field of 3-5 gauss, 
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Only the four transitions shown in fig. 1 were studied in detail, however, 
and the H.F.S. intervals and interaction constants which resulted after 
allowing for the field dependence were 

vy = 2884-7 + 0-2 Me/s, 

vg= 2171-5 + 0-1 Me/s. 
The interaction constants A and B with these intervals are 

A= — 446-97 + 0-04 Me/s, 

B= — 304-25 + 0-83 Me/s. 

Substitution of these intervals and the AF =0 frequencies v, and Ve 
into eqns. (1) and (2) gives 1-:6433+0-0002 as an accurate value for gy. 
The same result was obtained from both eqns. (1) and (2), thus showing that 
third-order perturbation theory yields a sufficiently accurate expression 
for the energy levels in this case. No direct measurement of vs was 
carried out. 


The observed frequencies v, and v, for the AF =0 transitions in different 
fields H,, and the energy intervals », and », calculated using 
eqns. (1) and (2) 


Observed frequencies Calculated intervals 


(39K) 
(Mc/s) B ] 
(Mc/s) (Mc/s) 


20-00 
30-00 
40-00 
60-00 
120-00 


§ 5. HyPERFINE STRUCTURE IN INTERMEDIATE COUPLING 


Atomic bismuth has three electrons missing from the nearest closed 
shell and the five states of lowest energy belong to the p® configuration, 
for which the terms allowed by the Pauli principle are *Psj., ?Pyj., 7Ds/, 
2D,/. and 4S,,. in pure LS coupling. The coupling is intermediate between 
LS and jj for J=3/2 states as indicated in fig. 4, so we represent the 
wavefunctions of these three states as linear combinations of the form 

ob(S = 3/2) =cyuy + Cotta + Cgg = dV, + dv, + dgvs 5 Ce Ce) 
where the v;, are the wavefunctions associated with the states k=1, 2 and 
3 in pure LS coupling and the wu, are the corresponding wave functions 
for pure jj coupling. The magnitudes of the c;, and the d;, depend on 
relative strengths for the Coulomb-repulsion and spin-orbit terms in 
the Hamiltonian, and there is a different set for each of the J =3/2 states. 


P.M. 4T 
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The wavefunctions for the J =1/2 and J=5/3 states are independent of 
coupling. The state of coupling is described by a parameter B, which can 
be determined by comparing the observed fine-structure with the theory, 
and the coefficients in eqn. (4) depend only on the energy W,,, for each state 
k, and B. Formulae for the calculation of the c's are given by Breit and 
Wills (1933), Inglis and Johnson (1931) deal with the d’s, and the general 
theory of intermediate coupling is discussed by Condon and Shortley 
1931). 
a the ground state of bismuth, one finds that X=3111-7cm™, 
B=3-10, c,= —0-188, cy= +0°345, cz= +0-918, d,= —0-354, d,= —0-515 
and d,;=+0:779. The extent to which which the parameters 8 and X 


Fig. 4 


(K=1) 


ENERGY PARAMETER G 


O | 2 3 4 
INTERMEDIATE COUPLING PARAMETER 8 ——> 


Atomic fine structure due to three p electrons in a state of intermediate coupling. 
The parameter £, defined by 


p2 _ 40,(w,+ 3)(w,—2) 


9(w,+3/5) : 
is a measure of the coupling, the value being 0 for pure LS coupling and 
oo for pure jj coupling. The ordinate w=2 a » & convenient 
j= 


function of energy, is plotted for the range B=0 to 4. The points 
represent optical results plotted at the value of the abscissa which 
affords the best fit to the theory. The states J=1 /2 and 5/2 have g,, 
A and B values which are independent of 6, although they may be 
perturbed by configuration interactions. 
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describe the energy levels is shown in fig. 4; the deviations may be due to 
the admixture of other configurations. 


5.1. Calculation of g, 


The operator g, = (g, Lz +9587) is diagonal in an LSJ M representation 
and the values of (LSJ |g,|LSJ) are well known, so it is convenient to use 
the LS wavefunctions u, for the calculation of g,. First order theory 
gives 

J. = Ay? (10419, [41 ) + da? (Wal j\U2) + ds? (uslgz|es >, \ (5) 
gj( J = 3/2) = 4,79, (?Psj2) + de°9,(?Dsj) + d3°g,(4Sz/0)- 


Substituting g,(?Dg/.)= — 0-799, gx Psi 2) = —1:334, g,(4Ss/.) = 2-002 and 
the d alten already listed into (5) gives g, = 1-651 for the bismuth ground 
state, in reasonable agreement with the experimental result — 1-6533(2). 


5.2. Calculation of A 


In order to calculate the expectation value of any operator V=>V (rj), 
r; being the position vector of the jth electron, one can use the jj wave- 
functions v;, reduced to the form of sums and products of the single-electron 
wavefunctions (/, 7, m), since the matrix elements <1,j,m|V|l',j’,m’) can 
be calculated using the Dirac equation. Combining results for the p?* 
configuration with (4) gives the magnitude of (V) for the state J =3/2, 
M =3/2 in terms of the single-electron matrix elements. 


CV) B= CRAIVIBE) + (ea? 4) DIV IRE) + (cat dex?) 8417182 
+ (£5? + c37) (44[V 54) + (2/-v/5)ca(ca — 2¢) 
(P= HV Rb) ( (2/4/5) C9(2cs — 61) (33|V [32 ) 
+ ($02? + cs") (3 —$|V [3 -4).- chro sO) 


The single-electron matrix elements of hA =2i(uy. H,), where H, is 
the magnetic field at the nucleus due to the jth electron, have been 
calculated by Breit and Wills (1931), who show that 


ha'= = + 29 ag, <r ‘avis l's ‘(1Z;) =h3(33|A[3$), 

ha" = + 297g, (1 dave F, AS hs<33|A[35 ), 

ha" = — 29°9 <7 daveG, =hy/ 234 [33 >, (7) 
A(SJ = $) = (ey? + $C" + C3 Na" 5s ie a" — 44/3 C9(Cy — C3)a", 


i 
A(J =3) =0-8a' + 0-2a”, 
A(J=4) =a", 
where F,’ (1Z;), F,”(1Z;) and G,(/Z;) are relativistic factors of the order unity 


which are tabulated conveniently by Kopfermann (1958), page 445. 
One way of estimating (r~*)ay is to make use of the expression 


G& =BX = 212Z,(r-8 av (IZ;) Segre (5) 


where @ is the fine-structure interaction constant, which depends on the 
energy separation between two states independent of coupling, and the 


472 
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coupling parameter f (fig. 4). The effective atomic number Z; is approxi- 
mately equal to Z—4 for p electrons. Taking X=3111-7 em; B= 3°10; 
the constants c already listed and the nuclear magnetic resonance value 
for g; we get +456 Mc/s for A(J=%). This is opposite in sign to the 
experimental result — 446-97 40-04 Mc/s. The eqns. (7) also predict a 
value wrong by a factor of five for the ratio a” /a’, which can be deduced 
from the optical measurements of A(J = $3) and A(J =>) made by Zeeman 
et al. (1930). ; 

Breit and Wills (1933) suggested that such discrepancies might be due 
to interactions with other electronic configurations. Similar inconsistencies 
in the ratios of parameters have been found in the case of the 4s4p 
configuration of gallium (Koster 1952) and the 5p configuration of rubidium 
(Phillips 1956) and these have been explained in terms of s electron 
excitation. Koster (1952) showed that to a first-order approximation 
the fine structure is not affected by the admixture of excited s states, 
though it may have a significant effect on the H.F.S., so it may still be 
possible to use eqn. (6) to estimate A(J= 3) if the a’s are treated as 
parameters to be found by experiment. Ifa’ and a” are deduced from the 
optical H.F.S. measurement of A(J=4) and A(J=3), thus making 
first-order correcting for the configuration effects, and a” is estimated 
from a”, the more accurately known parameter, A(J=3) comes to 
—413+90 Mc/s, where the error reflects the uncertainty in the optical 
measurements. This result is in reasonable agreement with the experi- 
mental value. Further evidence for such configuration interactions in 
similar atoms is given in the following paper on the antimony isotopes. 


5.3. Nuclear Quadrupole Moment 


In order to calculate @ from the measured quadrupole interaction 
constant B it is necessary to estimate q,, which is defined by 


dy = (— >; 3(3 cos? —1)),, : 

and then use hB=e?¢,Q. Schuler and Schmidt (1936) calculated the 
single-electron matrix elements of q, and substituted them into (6) to get 

qy(J = 2) = — BCP dav LR, (1Z;) (61? — 5?) + V/ 28, (1Z;)co(ey+¢3)]. (9) 
The relativity correction factors R(/Z,;) and S(/Z;) are listed in Kopfermann 
(1958, p. 448) and the c’s have already been given. Koster (1952) showed 
that the admixture of some excited s state does not affect the fine structure, 
and likewise, it should not contribute to the quadrupole interaction, so 
eqn. (9) may give a reasonable value for gq when (7-3 )ay is estimated from 
fine-structure data and eqn. (8). When this is done one finds Q = — 0-340. 
After making a Sternheimer (1952) correction, the result is increased 
slightly to —0-35b. The error in Q due to the uncertainty in Bis +0-001b, 
but a much larger error is expected in Q because of the effects of other 


configuration interactions on the estimation of q, and the uncertainty 
in the Sternheimer correction. 
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ABSTRACT 


The hyperfine structure of the J=3/2 ground state of 1248b and 1238b has 
been studied by the atomic beam magnetic resonance method, and the 
results are consistent with the following values for the interaction constants: 


9 = —1-9705 + 0-0002, Ay4= — (299-0344 0-004) Me/s, 
A y23= — (162-451 +.0-003) Mc/s, By; —(3+68 + 0-02) Me/s, 
and Bys3= — (4-67 + 0-03) Me/s. 


To the accuracy of the measurements there was no evidence for an octupole 
interaction. 

From the B values one obtains Q,.,;=—(0:20+0:03) x 10-24 em? and 
Q123= —(0-26+0-04) x 10-24 em?, and comparison of Aj,.,/A1.3 with the 
ratio of g42;/9123 already published gives —0°318+0-003% for the hyperfine 
structure anomaly. Intermediate coupling theory accounts satisfactorily for 
gz, but an admixture of excited s states is necessary to explain the observed 


magnetic interaction and anomaly. 


§ 1. INTRODUCTION 


LOWENTHAL (1929) investigated the optical Zeeman effect of the SbI fine 
structure and showed that Russell-Saunders coupling of the valence 
electrons was not complete. He assigned an experimental value to the 
ground state g;. Crawford and Bateson (1934) studied the optical H.F.S. 
of SbIV and assigned nuclear spin (/) and nuclear magnetic moment (17) 
values to the stable 14Sb and !3Sb isotopes. Loomis and Strandberg 
(1951) used the microwave resonance method to obtain an accurate 
value for the ratio of the quadrupole moments of these two isotopes 
and so stimulated further experiments on the optical hyperfine structure 
of SbII by Sprague and Tomboulian (1953). Murakawa (1955) using 
separated isotopes, and Gordy (1955), obtained results for the ratio of 
the quadrupole moments somewhat different from those of Sprague. 
No direct optical observation of the atomic ground state H.F.S. has yet 
been made. 

Dehmelt and Kriiger (1951), and Wang (1955), measured the ratio 
Q101/Q123 to about 1 in 6000 using pure nuclear quadrupole resonance 
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techniques, and the ratio j1491/4123 has been measured accurately in nuclear 
resonance experiments by Cohen et al. (1950), Collins (1950) and Proctor 
and Yu (1950). Hisinger and Feher (1958) have recently re-measured 
L491/H#193 and used the Electron Nuclear Double Resonance (ENDOR) 
technique to measure the ratio Aj.;/A493 for the simpler s,/.-like electronic 
configuration of Sb in silicon. 

We have used an atomic beam magnetic resonance apparatus to 
investigate the H.F.S. of the ground states of Sb and Sb. The 
field-dependent AF=0 transitions and all the AF=1 transitions were 
observed with both isotopes and from the results we deduce values of 
97, Q121 and Q4.3, and also the H.F.S. anomaly, which is defined by 


121 At93 = (@421/4193)(9193/9121) — 1. 


§2. APPARATUS 


The apparatus for this experiment was the same as that used in the 
study of 2°Bi by Title and Smith (1960) (see preceding paper, which 
will be referred to as I). 

Antimony dissociates at a higher temperature than bismuth so the 
oven snout was lengthened and made of tungsten instead of molybdenum. 
With this modification a convenient beam intensity was obtained when 
the end of the snout was hot enough to bring about a high degree of 
dissociation. Residual molecules were intercepted with a stop wire, and 
the signal-to-noise ratio was further improved by modulating the beam at 
14 ¢.p.s. instead of 25 c.p.s. The output from the phase-sensitive detector 
was fed into a pen recorder which traced the resonance as the oscillator 
was scanned through the resonant frequency. 


§ 3. PROCEDURE 
The energy level diagrams for 124Sb (J =3/2, [=5/2) and 12Sb (J =3/2, 
I =7/2) are shown in fig. 1. The levels are inverted because A is negative. 
The interaction constants are related to the zero field intervals as follows: 


LISb 123Sb 
vy=4A +4B/5 + 320/5 v3=5A+5B/7+40C0/7 
v= 3A — 9B/20—84C/5 ¥9=4A —2B/7—880/7 
v3=2A —4B/5 + 960/5 v3 =3A —5B/7 + 88C/7 


where C is a possible magnetic octupole interaction constant. 

v, and vz are the only two AF=0 transitions which satisfy the 
refocusing conditions of the machine. They were followed up to fields 
of 100 gauss, and the magnitudes of the constants g,;, A and B deduced 
therefrom were used to calculate the approximate frequencies of the 


AF=1 transitions. These were then measured directly to give more 
accurate values for A and B. 
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Since m,=0 levels are field independent to the first order, it is desirable 
to observe AF = 1 transitions involving these levels because small fluctua- 
tions in magnetic field during the experiment are less important, the 
field correction to the final result is smaller, and small field inhomogeneities 
in the C magnet do not deform the line shape. All the refocused AF =1 
transitions were scanned at C fields of about 1 gauss, and vp, the only field 
independent resonance which satisfied the refocusing conditions, was 
measured accurately. A double resonance technique introduced by Davis 
et al. (1949) was used to observe the field independent transitions v, and v3. 
Two triode cavity oscillators were coupled to the same r.f. loop, the first 
was set on one of the two maxima associated with the v, o-transition 
(Woodgate and Hellwarth 1956), while the second was scanned through 


Fig. 1 
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the frequencies expected for v, and vs. The refocused beam intensity 
was reduced on resonance because of transitions to the upper or lower 
level and the double-humped y, and vz resonances traced by the recorder 
thus appeared inverted. Careful adjustment of both power levels was 
essential. The procedure was then repeated with the first oscillator set 
on the other maximum of the v, resonance, but no shifts of frequency 
were found. From the three zero-field H.F.S. intervals measured in this 
way we deduced A and B and showed that C is zero within the accuracy 
of the measurements. The results are given in table 1. 
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The sign of A;),; was determined by a method used previously by King 
and Jaccarino (1954). Atoms with one sign of m,, say m,;= +4, were 
removed by a stop wire in the A magnet and the low frequency AF’ =0 
transition v, was used to refocus the m;= —+ atoms on to the detector. 
A second frequency v, (fig. 2) was then applied to a separate r.f. loop 
some two inches nearer the detector, and the refocused beam decreased 
in intensity, showing that A is negative; had A been positive there would 
have been no change in intensity. A similar effect was obtained when 
the frequencies vz, and y, were used in the reverse order and the atoms 
with m,=—4 were suppressed instead of those with m,= +3. 


Table 1. Zero field intervals and interaction constants 
for 1218b and 123Sb 


121Sb (Me/s) 238b (Mc/s) 


vg=595-12+0-01 vg=484-02+0-01 
v,=819-454+ 0-01 vy= 648-46 + 0-01 
v,=1199-08 + 0:01 v,=815-60+ 0-01 
A = — 299-034 + 0-004 A= — 162-451 + 0-003 
B= —3-68 + 0-02 B= —4-67+0-03 


g;= —1-9705 + 0-0002 


Fig. 2 
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An accurate value of g, can usually be obtained by substituting the 
HES. intervals in the perturbation theory equations for the AF=0 
transitions. However, the g,’s deduced from v, and vg differed by more 
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than the probable error when perturbation theory up to the third order 
was used, so the secular equation was solved numerically using a digital 
computer. The values of g, deduced from each transition were then 
identical over a wide range of C fields. 


§ 4. Discussion or RESULTS 
4.1. Fine Structure and g; 


Atomic antimony has two 5s electrons and three 5p electrons outside 
a closed sub-shell and the coupling of the p electrons is intermediate 
between LS and jj. The ground state wave function can be written in 
the form of eqn. (4) of I with 


c, = 0°361 + 0-002, Cy= — 0°689 + 0-004, Cz = — 0-628 + 0-006, 
d,= —0-188+0-005, d,= + 0-057 + 0-004, d, = 0-980 + 0-002, 
a = 3270-5 cm, B=0-93 + 0-02. 


This value of §, which is slightly different from that obtained by Inglis 
and Johnson (1931) using the same theory, was estimated by making a 
simplified least squares fit of the theory to the fine structure measured 
optically by Lowenthal (1929). The error in f gives an indication of the 
closeness of the fit of the intermediate-coupling theory to the fine structure 
and the quoted errors in the c’s and the d’s are those which arise only 
from the uncertainty in f. It will be noted that the state is almost 
entirely LS coupled, being 98% 4S). 
The theoretical estimate, neglecting relativistic effects, of 


gy = —1-971 + 0-004 


is obtained by substituting the above values of d; in eqn. (5) of I. This 
is in good agreement with the experimental g, = 1-9705 + 0-0002. 


4.2. Comparison of A with Theory 


The intermediate-coupling theory fits the five fine-structure levels with 
only two parameters 8 and d@, and the experimental agreement is quite 
good, so one might expect this theory to account for the observed 
hyperfine structure in a satisfactory manner. However, the calculation 
of A,., using this theory (eqn. (7)) gives Ay.;=—15+5 Mc/s, where 
the error is due to the uncertainty in 8. For this calculation we used the 
g, measured by Eisinger and Feher (1958), and the values @= 3041-6 cm™ 
and (r-3),,,=70-8 A’, which were obtained by substituting 6 =0-93 and 
X =3270-4cm~ in (8). The experimental value of A,,, is — 299-032 Mc/s. 
It is possible that the discrepancy is due to configuration interactions 
which do not alter either g, or the form of the fine structure. In addition, 
the ratio a'/a” deduced from (7) should be 0-17 for a pure p® configuration, 
whereas optical measurements of A(P,j,) and A(D5/.) by Badami (1932) 
yield 0-22. If we try to allow for the effect of any admixed contigurations 


1296 P. GC. B. Fernando et al. on the 


by treating the optical results a’ = 1068 + 350 Mc/s and a” = 4740 + 140 Mc/s 
as experimental parameters in the intermediate -coupling theory, a method 
which gave a reasonable fit in the case of Bi, we obtain 
A(J =5/2)= + 274 + 400 Me/s. 

The error here is associated with the optical measurements and the 
corresponding error due to the uncertainty in f is only 7 Mc/s. No 
conclusion about the validity of this procedure can be drawn because 
of the large error in the predicted value of A(J =3/2). 


4.3. Hyperfine Structure Anomaly 


The ratio 9451/9193 has been measured by Eisinger and Feher (1958) 
and this may be combined with our measurement of Aj.;/Aj.3 to give 
(—0:317+0-003)% for the H.F.S. anomaly j.;Aj.3. Hisinger and Feher 
have also measured 45,Aj93 for an atomic configuration which is essentially 
a pure s state and they obtained (— 0-352 + 0-005)%, which is only slightly 
larger than that measured for the intermediately-coupled p* configuration. 
They have considered several nuclear models and conclude that the model 
which gives the best result for j445;, Qj2; and A(s) is the configuration 
mixing model, for which A(s)= —0-365%. 

The Breit Rosenthal contribution to A is ~ 10-* % for the atomic s state 
and is much smaller for p states. We have estimated the Bohr—Weisskopf 
anomaly for the configuration mixing model, using (6), and conclude 
that A(p?) is probably smaller than —0-08°, if off-diagonal terms are 
neglected. The anomaly we have measured is typically that expected 
for an s state, and a small s state admixture could also account for the 
magnitude of A. 


4.4, Configuration Interactions 


The only configurations which can be mixed with the antimony ground 
state to produce a large hyperfine anomaly and almost no change of g, 
are those in which s electrons are excited to higher s states. Koster (1952) 
and Phillips (1956) accounted for inconsistent values of a’/a” in Ga and Rb, 
similar to those which we now find in Sb, by means of outer s electron 
excitation. Abragam (1955), Heine (1957), Sternheimer (1957), Cohen 
et al. (1959) and Blinder (1960) considered the exchange interactions 
between the aligned valence electrons and s electrons in different shells, 
and showed that different spin-up and spin-down electron densities at 
the nucleus are associated with each of the s shells. The contact term in 
the expression for A is then a sum of several terms, each with a sign and 
magnitude dependent on the overlap of the s and valency electron wave 
functions, and the resultant A is sensitive to the s wave functions used, 
No detailed calculation of this sort has been carried out for antimony. 

The breakdown of LS coupling contributes only —15 Me/s out of the 
total observed value of — 299 Me/s to A,.,. If we assume then that most 
of the magnetic interaction is due to s electron excitation, any s state 
admixture which satisfactorily accounts for A gives automatically an 
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H.F.S. anomaly of the right sign and magnitude, irrespective of the 
detailed form of the calculations. In addition, excitation of s electrons 
to higher s states does not alter g,, since the contribution from these 
electrons remains zero. 


Table 2. Values for Q11, Y123 and Qyo1/Q193 reported in the literature 


Method Gia @123 Qi21/@ 128 
el Oss oni" x 10-24 cm? 

Dehmelt and Kriiger (1951) — — 0:7843 
nuclear quadrupole + 0-0001 
resonance 

Loomis and Strandberg 
(1951), microwave spec- — 0-35 to — 0-45 to 0-791 
trum of Sb H,D — 0-50 — 0-62 

Sprague and Tomboulian si ky le — 
(1953), optical H.F.S. 

Murakawa (1955), optical —0-52+0-1 —0-67+0-1 0-73 
HEF'S. + 0-06 

Wang (1955), nuclear _ — 0-7844 
quadrupole resonance + 0-0001 
in crystals 

Gordy (1955), microwave —0-:8 —1-0 0-784 
absorption + 0-002 

As reported in this paper | —0-20+0-03 | —0-26+0-04 0-784 

+ 0-008 


4.5. Quadrupole Moments 


Sternheimer (1957) has published a general discussion of the effects of 
core excitation with particular reference to the calculation of quadrupole 
moments. He has shown that the corrections can be quite large even 
when fine or hyperfine structure data are used to calculate (r~*),,, and 
that they are sensitive to the wave functions used for the core and valence 
electrons. If the excitation of s electrons to higher s states is the pre- 
dominant effect, however, both the quadrupole interaction and the fine 
structure remain almost unchanged and we can use eqn. (9) of I to estimate 
Q. Quadrupole moments calculated in this way are compared with 
results obtained by other methods in table 2. The errors assigned to our 
Q values are those associated with the uncertainty in f, and no Sternheimer 
corrections have been made. Although there is some considerable dis- 
agreement about the magnitudes of the quadrupole moments, the signs 
are all negative and our value for Qj49;/Qie3 agrees with the quadrupole 
resonance results. 
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On the Scattering of Phonons by Spins at Low Temperatures 
Experimental 


By H. M. Roszenpere and B. Susaxt 
The Clarendon Laboratory, Oxford 


[Received August 29, 1960] 


THIS paper reports some preliminary experiments which show how 
the heat conductivity of a dielectric crystal may be profoundly influenced 
by the presence of a small number of paramagnetic ions. Measurements 
of the thermal conductivity between 2 and 20°K have been made on 
crystals of Analar grade zine sulphate, ZnSO,.7H,O (maximum impurity 
0-006%) and also on zine sulphate crystals which contained 0-25, 1 and 
5-6at.% of FeSO,.7H,O. These measurements were compared with 
those on crystals of zine sulphate which contained 5-7at.% of a non- 
paramagnetic impurity, MgSO,.7H,O. All the crystals were 2 to 3cm 
long and had a square cross section with a side of 3to4mm. The measure- 
ments were made in a manner similar to that described by Rosenberg 
(1955). 

Some of the results are shown in the figure (a). The curves for the pure 
material are on the whole of the usual type which one expects for the 
thermal conductivity of a dielectric crystal (see Berman 1953, Klemens 
1956, or Ziman 1960, for general reviews). The rapid rise in conductivity 
on the high temperature side of the maximum is due to the decrease in the 
effectiveness of the scattering of the phonons by impurities or isotopes, 
and the decrease in conductivity at temperatures below the maximum is 
caused by the scattering of the phonons by the boundaries of the specimen 
or of small crystallites. It is unfortunate that it is not easy to grow very 
perfect crystals from a saturated aqueous solution and there is considerable 
scatter (from 25 to 50%) in the data from different crystals of the same 
composition. 

The main interest lies in the measurements taken above the conductivity 
maximum. The extra thermal resistance in the specimens containing 
Fe?+ is shown in the figure(b). It will be noted that the addition of 
5-6% Fe?+ decreases the conductivity and completely changes the 
temperature dependence from 7’ to about 7~°°. This change cannot 
be ascribed merely to an impurity scattering effect, because as the 
figure (a) shows, the conductivity of a specimen containing 5-7% Mg is 
very much closer to that of the pure material. This can be emphasized 
still further when we note that the thermal resistance due to impurity 
scattering is proportional to the square of the mass difference between the 


+ On leave from Zaklad Niskich Temperatur Instytutu Fizyki Polskie), 
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impurity and the other atoms (see Ziman 1960). Since the mass difference 
between Mg and Zn is about 41, whereas between Fe and Zn it is only 9:5, 
the effect of the Mg impurity should be very much greater than the Fe 
impurity instead of being much less. 
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Discussion 


We suggest that the large decrease in conductivity which is observed 
for the specimens containing 5-6% Fe?+ is due to Raman type interactions 
between the phonons and the low-lying spin levels of the ferrous ion By 
this we mean that a phonon interacts with a spin, flips it to another level 
and another phonon is emitted. Since, for the Fe?+ ion, the spin-lattice 
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relaxation time is thought to be very short (see, for example, Tinkham 
1956), the probability for this mechanism will be high. The theory of 
these processes has been given by Orbach (1960) and is presented in the 
following letter. He has also calculated the scattering due te direct 
processes, i.e. those in which a phonon of just the right energy is absorbed 
so that a spin may be flipped without the need for the emission of another 
phonon to maintain the energy balance, but he finds that the probability 
of such processes is much lower in the temperature region of our experi- 
ments than are the Raman processes. 

Unfortunately a detailed calculation of the effect is not possible because 
the splitting of the lowest levels of the Fe?+ ion is not known. Some 
specific heat measurements have been made on the 5:6% Fe?* crystal, 
however, and these have an anomaly with a peak at about 1-:8°K, which 
can be ascribed to a splitting of 4:3+40-2 degrees. If such a splitting is 
assumed, then Orbach’s calculations yield an additional thermal resistivity 
of about the order of magnitude which we observe and with a temperature 
dependence which seems to vary with concentration as in fig. 1(b) (see 
figure in following letter). 

It is clear that moré experiments must be made using other types of 
crystal (preferably without water of crystallization, because these are 
usually much more perfect). It is just possible that the difference in the 
curves which we observed are due to the fact that some specimens are not 
very perfect. It is known that some hexahydrate can co-exist with the 
heptahydrate and that there are in fact two crystallographic forms for the 
heptahydrate (Giauque et al. 1950). It does not seem likely, however, 
that the large change observed for the specimens containing 5-6% Fe?+ 
can be due entirely to defective crystals. Other types of spin-phonon 
interaction in heat conduction experiments have in fact been observed in 
ferrites (Friedberg and Douthett, 1958) and antiferromagnetics (Slack 
and Newman 1958) and these give further grounds for confidence in our 
results. 
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On the Scattering of Phonons by Spins at Low Temperatures 
Theoretical 


By R. OrBacu 


University of California, Department of Physics, Berkeley, California 
and 
Clarendon Laboratory, Oxford 


THE previous paper (Rosenberg and Sujak 1960) has described the observed 
effects of paramagnetic impurities on the thermal conductivity of a 
dielectric salt: FeSO,.7H,O in ZnSO,.7H,O. It is the purpose of this 
letter to briefly analyse their results in terms of a theoretical model. 

In the usual treatment of spin-lattice relaxation in paramagnetic salts 
(Van Vleck 1940) there are two important relaxation processes which are 
considered. The first, the so-called ‘direct’ process, involves the emission 
or absorption of a single phonon of the energy corresponding to the splitting 
kA (whether due to Zeeman or crystalline field effects) between the spin 
states of the ion. For Zeeman splittings the phonon energy is small in 
moderate fields, ~0-1°k for a g=2 and a field of 1000 gauss. There is, 
however, a second relaxation process, the so-called ‘Raman’ process, 
wherein a phonon is inelastically scattered by the spin and it is the difference 
between the initial and final state phonon energies which is equal to the 
splitting of the spin states KA. The latter process involves all phonons, 
but primarily those at energy ~kZ' because they are most abundant. 
If kT> the splitting, then the Raman process is usually dominant in 
spin-lattice relaxation. 

We now reverse the tables and ask how the thermal resistance of a 
dielectric salt will be affected by paramagnetic impurities. That is, we 
must calculate the scattering cross section for phonons in the presence of 
spin—phonon interactions. We can then find the expected thermal 
resistance. 

It is evident that if the splittings of the spin states A<7, the direct 
process will lead to a negligible thermal resistance as the great bulk of 
phonons which carry the heat will have energy ~k7' and could not be 
scattered by processes which are limited to phonons of the splitting energy. 
Thus, thermal resistance due to direct processes will be important for small 
splittings only at very low temperatures, on the order of <A. However, 
as Raman process scattering affects all phonons, it can be expected to be 
the dominant mechanism for providing a significant thermal resistance at 
higher temperatures. 

To gain a semi-quantitative estimate of the effects of Raman processes, 
we construct a phenomenological Hamiltonian # coupling the spin 
system to the lattice (Orbach 1960), 


is = iB >». j, mi SS? + S/S?) Eom mje? e ° . ° (1) 
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where [ ig a phenomenological coupling constant, S/ is the th component 
of the spin at the [th lattice site measured in units of 4, we sum over all 
spin sites J, and «, is i, jth component of strain for the phonon (a). We 
restrict j to equal either a or y in order that #’ can ‘flip’ aspin. Equation 
(1) is then an expression for a Hamiltonian which annihilates or creates 
phonon (1), flips a spin, and creates or annihilates phonon (2). It applies 
to spin states which are not Kramer's degenerate in that it predicts a 
seventh power temperature dependence for the relaxation time. 

If we consider Fe2+ in ZnSO,.7H,0O, the Fe? will be in a predominantly 
tetragonal crystalline field with a slight orthorhombic distortion. Together 
with the spin-orbit coupling, this results in the removal of all orbital and 
spin degeneracy as the number of electrons is even. Normally, with only 
a tetragonal field, one would expect a doublet lowest, but the presence of 
the small orthorhombic term will slightly split the doublet levels. The 
specific heat measurements of Bailey (unpublished) show that under 
these considerations, the splitting of the lowest levels (presumably here 
the doublet levels) is Ax4:3°K. Such a level scheme implies that eqn. (1) 
has the proper form. 

We compute the spin-lattice relaxation time 7’, from eqn. (1) by standard 
means (Orbach, loc. cit.) 


1/7, = 908412 /(ar3p2v!) x (kT hy. dike aborted 


Here, p is the density, v the phonon velocity (assumed the same for all 
modes for simplicity), and S the spin value. Likewise, one can easily 
compute the relaxation time for phonons via the interaction Hamiltonian 
(1) (Orbach, loc. cit.), 


1/r = S4P2N ,/(12mp2v7) x {w2[w? + (kA/h)2} =. . (8) 


where there are V, spin impurities per unit volume, A is the splitting 
between the spin states in °K and w is the phonon frequency. From a 
comparison of (3) and (2), it is clear that a short spin—lattice relaxation 
time 7’, implies a short phonon mean free path= vr. 

If we compute the thermal resistance due to eqn. (3) by the usual methods 
(Klemens 1958) neglecting questions of thermal equilibrium within the 
phonon spectrum, we find a thermal resistance W given by, 


1/W = 408kp?|(mS4T2N ,) x (i/ kT)? 
x | wo? dw exp (he/kT)/{{co® + (kA/h)?] [exp (alk) 12}. (4) 
For A<T, this reduces to, 
1/W ~ 4p%kv®/(mS4T2N,) x (kA). =... . (5) 


From paramagnetic resonance line broadening due to spin-lattice inter- 
actions for Fe*+ in Zn, 'Tinkham (1956) found that T,~5 x 10-*see at 
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T'=20°K. We correct for the difference in sound velocities between ZnF, 
and ZnSO,.7H,O and, though the crystalline field effects are quite 
different, for an order of magnitude estimate we assume the spin-lattice 
coupling constant TI‘ stays roughly the same. Then, from (3), we estimate 
[?, using p=2g/em*, v=2-1 x 10°cm/sec, and an ‘effective’ spin S=#. 
We find I?~10-ergs?. Using this value, the thermal resistance W 
from eqn. (5) with, say, 5-6% Fe?+ in ZnSO,.7H,0 [i.e. N,=2-93 x 1020 
spins cm~*] 30 watts~! cm deg, in very good order of magnitude agree- 
ment with the mean experimental value (see figure). According to eqn. 
(5), W should be temperature independent and directly proportional to 
the splitting A and the concentration JN. 

In fact, for 5-6% Fe?+ in ZnSO,.7H,O, Rosenberg and Sujak find a 
thermal resistance which is temperature dependent (see figure). Further, 
the resistance does not appear to be directly proportional to concentration 
and the temperature dependence seems to vary with concentration. Both 
these apparent discrepancies can be removed by noting that our calculation 
of W may have been too naive. If we let A+0, then eqn. (3) will have 
the same frequency dependence as one would find from impurity scattering, 
and the thermal resistance as given by eqn. (5) would vanish due to the 
heat flux of those phonons whose energy is <7’ not being significantly 
attenuated. In fact, as A is significantly less than 7’ in the experiments 
of Rosenberg and Sujak, it would appear that just such a spurious effect 
might be entering here. Two apparently equivalent treatments have 
been given (Callaway 1959, Berman et al. 1959) to remedy such a situation. 
Both consider that the usually neglected momentum conserving (J) 
processes tend to bring the low frequency phonons back into equilibrium 
with the higher frequency phonons (of energy >k7') whose heat flux is 
much more effectively impeded by the non-momentum conserving 
processes (in their case, impurity scattering). We follow the treatment 
of Callaway because of its simplicity. We assume that a phonon of 
frequency w has a lifetime ry due to N-processes given by 


Lye Aa) Saga eh, Pen e+ <1 (8) 


where A isaconstant. This is the usual expression for transverse phonons. 
Longitudinal phonons would obey a law like 1/ry%Aw°(kT/h)? for an 
orthorhombic crystal, but, as we are interested in those phonons for 
which w<kT'/h, the N-process transverse phonon relaxation time (6) is 
much shorter than the N-process longitudinal phonon relaxation time. 
We make the rather reasonable assumption that the longitudinal phonons 
are in equilibrium with the transverse phonons and hence consider (6) as 
the overall effective N-process relaxation time. Our expression for 
the thermal resistance W is then, using the formalism due to Callaway, 


L/W =40%kp?/(nS4T2N,) x (A/T? 
: | w8 dw exp (he/kT)|{[o? + (b/w + a 7/A)(KL'/A)*] [exp (Hieo/kT) — 1} 
Be ei) 
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where « is a dimensionless constant giving the relative strength of the 
momentum conserving processes to the non-momentum conserving 
Raman processes. Explicitly, 

a= A(kA/h) x 127p*0' | (S*1T AN). owes es, 
Using the approximate theory due to Orbach (loc. cit.), we find that 
A~10-*4sec* utilizing the appropriate elastic properties of ZnSO,.7H,O. 
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The calculated thermal resistance W due to substitutions of 0-25, 1-0 and 
56 at. % ferrous sulphate in zine sulphate given by eqn. (7), plotted 
against temperature. The crosses are the experimental points of 
Rosenberg and Sujak (previous paper) for the corresponding salts. 


We evaluate the expresion for « for this value of A and for the values of 
p, v, A and I? given above. Respectively for 5-6%, 1% and 0-25% Fe2+ 
in ZnSO,.7H,O, a=0-001, 0:0056 and 0-022. Inserting these values 
in eqn. (7) and evaluating the integral, we find the expected thermal 
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resistance due to spin—phonon interactions in this salt. If we divide our 
order of magnitude results by a factor of ~ 2 in order to fit the experimental 
magnitude of W for 5-6% FeSO, at 10°K we find the theoretical curves of 
W against temperature 7’ shown in the figure. These are to be compared 
with the experimental results of Rosenberg and Sujak given in the pre- 
ceding paper which are also shown in the figure for convenience. Without 
any further adjustments the temperature dependence of W as predicted 
by eqn. (7) seems to fit the experimental data rather well for the sub- 
stitution of 5-6% Fe?+. The predicted change in temperature dependence 
with concentration and the relative magnitude of the thermal resistance 
also seem to be in excellent agreement with the experimental results for 
1-0 and 0-25% Fe*+ over the measured temperature range. 

This agreement would seem to justify our interpretation of the experi- 
mental results. In particular, the choice of the phenomenological spin- 
phonon coupling Hamiltonian (eqn. (1)) to give the corresponding phonon 
relaxation time (eqn. (3)), and the use of the V-processes to restore thermal 
equilibrium amongst the phonon modes would both appear to be consistent 
with the experimental results. Further investigations are planned in 
other better behaved substances (see preceding paper). It is hoped that 
the values of [ and A can then be pinned down more closely and more 
accurate experimental results can be obtained. 
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The Spin of '74Sb and the Spin and Moments of !??Sb 


By P. C. B. Fernanpot, G. K. Rocuzsrer and K. F. Surre 


Cavendish Laboratory, Cambridge 
[Received August 19, 1960] 


THE spin of 14Sb and the spin and hyperfine structure of 122Sb have 
been measured by the atomic beam magnetic resonance method. We 
find I,,,=3, which verifies Steffen’s (1960) 8 decay scheme for 124Sb, 
and J,..=2, in agreement with Pipkin’s (1958) independent measurement 
made by the paramagnetic resonance technique. The Af’ =1 transitions 
5/2, —1/2<—>3/2, 1/2 and 7/2, —3/2<>5/2, —3/2 were among those 
observed in 17°Sb, and from these it was deduced that 


Ajog= +212-01+40-06 Mc/s and Byyg= + 8-7 40-5 Me/s. 


The assignment of the positive sign to A,., was made on the basis of the 
determination of the sign of ,,,, by Hisinger and Feher (1958), that of 
Hig by Pipkin (1958), and that of A,., by Rochester (1960). The ratio 
Aj49/A1o, is consistent with the ratio 9495/91., derived from Pipkin’s 
measurements, but we cannot deduce a value for the hyperfine structure 
anomaly because the error of +1% in this g ratio is much greater than 
a typical s state anomaly (Fernando 1960). The uncorrected quadrupole 
moment deduced from Bj. is Qj29= (+ 0°47 + 0-03) x 10-74 cm?. 
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Activation of Single Phonon Infra-red Lattice Absorption 
in Neutron Irradiated Diamond 


By S. D. Surrx and J. R. Harpyt 
J. J. Thomson Physical Laboratory, University of Reading 


[Received October 13, 1960] 


Lax and Burstein (1955) have proved that the symmetry of the perfect 
diamond structure forbids single phonon absorption of radiation. They 
have also shown that such absorption may be allowed when imperfections 
are present, since these can provide a local coupling between the 
fundamental lattice frequencies and the radiation field. They have 
ascribed the specimen dependent infra-red absorption, shown by natural 
type I diamonds, to a coupling of this kind. Such an absorption spectrum 
is shown in the figure (a) and one can see that the continuous absorption 
commences near the maximum lattice frequency which is known from 
Raman scattering measurements. 

We have succeeded in activating a corresponding band in the 
absorption spectrum of a type II a diamond by fast neutron irradiation. 
Absorption spectra taken before and after a fast neutron dose of 
3x 101’/em? in a hollow uranium slug of the Harwell BEPO pile are 
shown in the figure (b). That taken prior to irradiation shows the 
negligible absorption characteristic of type Ila crystals. After the 
irradiation the spectrum shows strong continuous absorption, commencing 
sharply at the Raman energy and extending to lower frequencies, which 
can only be explained satisfactorily as single phonon lattice absorption. 
This effect has not been observed previously, although a number of 
workers (see Kaiser and Bond 1959 for references) have examined 
specimens known to contain appreciable numbers of radiation produced 
defects. Moreover it seems that it is relatively easy to activate this 
absorption in diamond, compared with silicon, for which Fan and 
Ramdas (1959) have reported a dose of 10!® fast neutrons/cm? to be 
necessary to produce a similar effect. This conclusion is confirmed by 
the fact that a control specimen of silicon irradiated with the diamond 
showed no measurable single phonon absorption afterwards. 

Kaiser and Bond (1959) have shown that the strengths of the 0-159, 
0-148 and 0-139 ev peaks in the type I spectrum are directly proportional 
to nitrogen content. They have suggested that these peaks are due to 
absorption by molecular stretching vibrations of the C-N unit. If one 
compares the two spectra in the figure it is evident that the only marked 
difference, within the regions of continuous absorption, is in the vicinity 
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of the two strong peaks (X) of the type I spectrum. We see that the 
0-148 and 0-139 ev peaks are present in both spectra, although the 
position of the second peak is slightly influenced by the strength of the 
first. They are not dependent on the presence of nitrogen and therefore 
cannot arise from vibrations of C-N units. Their presence and that of 
the 0-125 ev peak in spectra activated by two very different coupling 
mechanisms implies that all three arise from maxima in the density of 
plane wave normal modes of the host lattice, since this is the only common 


40 
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Absorption coefficient versus energy for: 
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The various maxima are assigned to those in the density of normal 
modes thus : T.O. denotes that arising from the maximum density of 
transverse optical modes, L.O., that due to maximum density of 
longitudinal optical modes, and L.A. that associated with the maximum 
density of longitudinal acoustic modes. 
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factor involved. However, the 0-159 ev peak is apparently characteristic 
of the nitrogen impurity. Even so, one cannot regard it as arising from 
a local C-N vibration since the frequencies involved lie within the band 
of allowed lattice frequencies and energy absorbed near nitrogen 
impurities must propagate throughout the crystal, and this is exactly 
the process envisaged by Lax and Burstein (1955). Either this peak 
is due to exceptionally strong coupling to relatively few modes, or it 
corresponds to a sharp maximum in the density of normal modes of the 
whole lattice. We believe that the second explanation is correct, since 
one expects the transverse optical (T.O.) modes to show very little 
dispersion and to produce such a maximum. One way to explain the 
absence of this peak from the damage induced absorption spectrum is 
to assume that the coupling involved falls off rapidly close to the Raman 
frequency. The figure also shows isolated peaks lying above the maximum 
frequency for which lattice waves can propagate. However, Montroll and 
Potts (1955) have shown theoretically that vibrational modes localized 
about lattice imperfections would, if present, occurin this region. The dipole 
moment associated with such modes could cause the observed absorption. 
The peaks present in the spectrum of the irradiated diamond have 
evidently been produced directly by the bombardment and are presumably 
associated with specific defects (e.g. vacancies and interstitials). This 
appears to be the first time that such an effect has been observed in any 
material. 

The continuous absorption produced by damage appears to be too 
strong to be associated solely with isolated point defects. A lower limit 
to the number required is given by assuming that they provide a coupling 
as efficient as that due to nitrogen atoms in type I crystalst. If one 
then considers the 0-148 ev peak, common to both spectra, one finds 
at least 102° defects/em? to be necessary to produce the observed intensity. 
This is at least an order of magnitude greater than the number one infers 
from the dose. Hence we suggest that most of the damage induced 
activation arises from regions of local disorder which probably involve 
more than 10? atoms/cm?. 

Since we have associated the absorption peaks with maxima in the 
density of normal modes we have classified them in the figure according 
to the branches of the spectrum believed to be responsible (e.g. the 
0-148 ev peak is believed to arise from the maximum density of 
longitudinal optical (L.O.) modes). The 0-139 ev peak is also associated 
with the L.O. branch. It could correspond to the maximum density of 
longitudinal acoustic (L.A.) modes, but this implies that the 0-125 ev 


+R. J. Elliott (Clarendon Laboratory Dept., Ref. 23/60) has shown that, 
however the nitrogen in type I’s is present, its extra ‘electron is probably 
equally shared with the host lattice. The strength of the resultant dipole 
will then account for that of the 0-159ev peak. It seems unlikely that 
irradiation produced point defects provide any more efficient coupling 
mechanism. 
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peak is due to transverse acoustic (T.A.) modes. Since Phillips (1959) 
has concluded that the T.A. mode energy at the centre of a (111) zone 
face is 0-063 ev it seems unlikely that there will be a large number of such 
modes at twice this energy. Consequently if our assignment is correct, 
the 0-131 ev absorption minimum is associated with a similar feature in 
the density of L modes. 

Our suggestion that practically no coupling to the transverse optical 
(T.0.) modes results from damage appears to be in marked contrast 
with the silicon results of Fan and Ramdas (1959) since their strong 
20-5 micron band would seem, from the inelastic neutron scattering 
data of Brockhouse (1959), to arise from strong coupling to T.O. modes. 
We hope to clarify this by studying the activated absorption as a function 
of dose and annealing time. 

We may conclude that activation of absorption in the ‘ one phonon’ 
region by damage, strongly supports the theory that it is due to defect 
induced coupling to the fundamental lattice modes. Furthermore the 
common features of this and type I absorption indicate that the latter 
has a similar origin. Differences in detailed shape should thus be 
explicable in terms of the actual coupling mechanisms involved. 
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REVIEWS OF BOOKS 


Figures of Equilibrium of Celestial Bodies. By ZDENEK Koran. (University of 
Wisconsin Press, 1960.) [Pp. 135.] $3. 


INVESTIGATIONS into the shapes, gravitational fields, moments of inertia, oscilla- 
tions, etc. of approximately spherical masses like the planets have led to 
important work by Clairaut, Laplace, Stokes and many other experts. This 
branch of applied mathematics has come into prominence recently on account 
of the use of gravity surveying and satellite orbits in the study of the shape 
of the Earth. 

In the first of the three principal sections of his book, Professor Kopal 
presents a treatment of the problem on classical lines, based on the idea of a 
general disturbing potential, and using expansions in harmonic functions. In 
particular, the disturbing potential is taken as one representing the rotation 
of the body, or as one corresponding with equilibrium-tidal influences of an 
external mass. 

Professor Kopal then extends the treatment to the second order of the 
distortions. In the third principal section, he considers ‘ interaction pheno- 
mena’ connected with the essential non-linearity of the equations for the 
rotating body under tidal forces: at this stage, needless to say, the formulae 
become extremely complicated. 

Much of the advanced work set out in the book is the result of the investi- 
gations of Professor Kopal himself. 

There is also a short section dealing with dynamical considerations under the 
title ‘ non-radial oscillations ’. 

Formulae connected with the external potential of the Earth are of course 
fundamental in the geodetic applications of artificial satellites, but the book 
contains no treatment of satellite orbits, and the second part of the title referring 
to ‘ motion of artificial satellites ’ may therefore be a little misleading. ae 
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Fig. 1 
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Surface pits at the grain boundaries of a thin aluminium foil. x 20000. 
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Surface pits at the grain boundaries of a thin aluminium foil. x60 000. 


Fig. 3 Fig. 4 


Grain boundary of aluminium foil 


The boundary of fig. 3 after 
before ageing. x 40 000. 


ageing. x40 000. 
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Pits at both surfaces of a thin aluminium foil. x80 000. 
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Etch pits in a thin aged aluminium foil. x 4000. 
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Large etch pits containing smaller hexagonal 
pits on the basal surface of ice grown ice. 
from the melt. (Replica.) 


Replica of basal plane showing spiral steps A line of pyramidal etch pits along a 
on etch pits. (Replica.) low-angle grain boundary. (Replica.) 
Fig. 5 Fig. 6 


Etch channels on basal surface initiated 
by an included dust particle and 
terminating in etch pits. (Replica.) 


Pyramidal hillocks showing concentric 
and spiral steps on the basal face of 


an ice crystal grown from the vapour. 
(Replica.) 
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Elongated etch pits on the prism face of A high concentration of hexagonal etch 
an ice crystal grown from the pits on the basal surface of a plate- 
vapour. (Replica.) like crystal grown from the vapour. 

Fig. 9 Fig. 10 


Rectangular etch pits on the prism faces Etch channels radiating from an impurity 
of an ice needle. speck in a plate-like crystal grown 
from the vapour. Each channel con- 
sists of a double row of small etch 
pits which probably represent the 
intersection of dislocation loops with 
the surface. (Replica.) 


Fig. 12 


‘Rosettes’? produced by dropping lead 


Etch channels oriented along the <1120> 
shot on a basal surface of ice at 


directions on the basal plane pro- sal 
duced by the movement of disloca- —20°c. (Replica.) 
tions under mechanical vibration at 

—20°c  (Replica.) 
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Fig. 1 


Cross section of a 0-001 in. copper foil after an anneal of 2 hours at 800°c 
to show the formation of helium bubbles near the grain boundaries. 
x 550. 


(0) 
Two films of copper (a) and (b) through which 1-4 x 101 alpha-particles em-2 
(32 «amp hr.) have passed. 
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‘hot-stage ’ of the electron microscope. 
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A thin film of copper containing helium atoms after bombardment 
with 1-4 x 10” alpha-particles cm-?. 


Fig. 5 
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Fig. 5 (continued) 


grain -—_ 
boundary 


(e) 
r films similar to that in fig. 4 but after (a) a 2-hour anneal at 350°C 
Os in vacuo, (b) 2 hours at 450°c, (¢) L hour at 750°c, (d) as for (c) but 
with the dislocations out of contrast, (¢) 2 hours at 750°c, near a grain 
boundary (7 x 10'6 alpha-particles em~?). 
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Typical micrographs of the alloys (x 250). 
(i) 1-91 at. % chromium (F.R.1./B.M.1.); 
ii) 12-lat.% chromium (J.M./J.M.); 

(ili) 38-5 at. % chromium (F.R.I./U.S.B.M.). 
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